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Yang-Mills theories in four space-time dimensions possess a hidden symmetry which 
does not exhibit itself as a symmetry of classical Lagrangians, but is only revealed on the 
quantum level. It turns out that the effective Yang— Mills dynamics in several important 
limits is described by completely integrable systems that prove to be related to the 
celebrated Heisenberg spin chain and its generalizations. In this review, we explain 
the general phenomenon of complete integrability and its realization in several different 
situations. As a prime example, we consider in some detail the scale dependence of 
composite (Wilson) operators in QCD and super- Yang— Mills (SYM) theories. High- 
energy (Regge) behavior of scattering amplitudes in QCD is also discussed and provides 
one with another realization of the same phenomenon that differs, however, from the 
first example in essential details. As a third example, we address the low-energy effective 
action in a A/" = 2 SYM theory which, contrary to the previous two cases, corresponds 
to a classical integrable model. Finally, we include a short overview of recent attempts 
to use gauge / string duality in order to relate integrability of Yang-Mills dynamics with 
the hidden symmetry of a string theory on a curved background. 
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We had the privilege and pleasure of knowing Ian Kogan as a close friend. 

In our informal chats it was clear that, whatever he was doing, there always 
was a single central problem that kept his mind constantly busy - the prob- 
lem of confinement. All other projects he pursued - ranging from the Quan- 
tum Hall Effect to strings and quantum gravity - he considered from this 
general perspective. He viewed the issue of integrability as a unifying concept 
which could potentially explain many unsolved problems of strongly coupled 
theories and was very enthusiastic about this new development concerning 
four-dimensional gauge theories. If it were not for his untimely death, he 
would definitely have had a profound impact on the whole field with his 
usual passion for new concepts and trends. 

1. Introduction 

Understanding the dynamics of gauge theories at strong couphng is one of 
the outstanding problems in quantum field theory. It has been anticipated, 
since the classic work by Wilson [1] , that the Yang-Mills dynamics at strong 
coupling can be reformulated as an effective theory of Faraday lines (gauge 
field flux) which, in its turn, can be reinterpreted as a string-like theory. A 
natural framework for discussing the stringy representation is offered by the 
multicolor limit. In this limit, one can express observables as a sum over 
surfaces (strong coupling expansion a la Wilson) but explicit realization of 
this program in gauge theories is still lacking due to absence of an operational 
gauge/string correspondence. It is expected that in the multicolor limit, 
the path integral over all gauge field configurations will be dominated by 
a saddle point - the so-called "master field". Polyakov has suggested [2] 
that the master field can be found as a solution to a two-dimensional sigma 
model defined on a nontrivial background (Lobachevski space-time). He also 
demonstrated that such models possess an infinite set of conserved currents 
and, therefore, can be solved exactly. It has long been speculated that 
similar nontrivial algebraic structures and, in particular, integrability can be 
revealed in four-dimensional quantum gauge theories as well. In the last few 
years several plausible indications that this actually occurs have appeared 
and several intriguing connections have been established between quantum 
gauge field theories and integrable lattice spin models. Recent progress has 
been in several directions, which we describe here. 

At first glance, the connection between four-dimensional gauge theories 
and integrable systems may seem surprising. Indeed, in the latter case one 
is dealing with quantum-mechanical systems, which have a finite number of 
degrees of freedom and the same number of conserved charges. For such 
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systems the degrees of the freedom are just phase-space variables involved in 
the Hamiltonian and the "evolution time" has a literal physical meaning. In 
contrast, Yang-Mills theories in four dimensions are complex systems with 
infinitely many degrees of freedom which are not integrable per se. Complete 
integrability emerges as a unique feature of effective Yang-Mills dynamics 
in various limits. The relevant degrees of freedom, Hamiltonians and "evo- 
lution times" within the Yang-Mills theory are different in different limits 
and their identification is not a priori obvious. In this work we will elabo- 
rate on three examples which illustrate the general phenomenon of hidden 
integrability. They are 

• Scale dependence of composite (Wilson) operators in QCD and 
super- Yang-Mills (SYM) theory 

• High-energy (Regge) behavior of scattering amplitudes in QCD 

• Low-energy effective action in M = 2 SYM theory 

In all three cases, the dynamics described by an integrable model corresponds 
to the scale dependence of corresponding observables with the "evolution 
time" being identified as the logarithm of the relevant energy scale. In the 
first two cases one is dealing with quantum integrable models, while the 
M = 2 low-energy action corresponds to a classical model. The integrable 
systems that emerge in this context turn out to be related to the celebrated 
Heisenberg spin magnet and its offspring. This model was introduced by 
Heisenberg in 1926 and describes a one-dimensional chain of atoms with an 
exchange interaction. 



where S„ = {S^, Sn, S^) is the spin— 1/2 operator of the re— th atom in the 
chain of length L and periodic boundary conditions are implied: Sl+i = Si. 
The model (jl.lj) is completely integrable and its eigenspectrum was found 
in 1931 by Bethe by an innovative method which is widely known nowadays 
as the Bethe Ansatz. Much later, it was understood [3,4] that the original 
Heisenberg model (|l.lj) can be generalized to arbitrary spins while preserving 
complete integrability. The Hamiltonian of a completely integrable lattice 
model describing a chain of interacting spin— s operators was found to be [4] 




n=l 



(1.1) 



L 
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Here the operator Jn,n+i is related to the sum of two spins in the neighboring 
sites, = s(s + 1), and H{x) is the following harmonic sum 

2s-l ^ 

H(.x) = Y,Y^ = i^i'^s + l)-ij{x + l), (1.3) 

l=x 

where ^{x) = dlnT(x)/dx is the Euler function. For s = 1/2, the two- 
particle spin takes the values Jn,n+i = and Jn,n+i = 1- In that case, 
i?(0) = 1 and H{1) = so that the two-particle Hamiltonian H{Jn^n+i) is 
given by a projector onto J„,„+i = subspace, H{Jn,n+i) = 3 - S„ • S„+i, 
in agreement with Hl.l|) . 

Approximately at the same time as the model ()1.2p was formulated, QCD 
calculations of the anomalous dimensions of twist-two Wilson operators [5] 
and high-energy asymptotics of scattering amplitudes [6] led to expressions 
involving the very same combination [V'(</) — "0(1)] with J being the Lorentz 
and conformal SL[2) spin respectively. In both situations, the appearance of 
the -i/^-functions is a generic feature related to the existence of massless vector 
fields (gluons). For almost two decades, this similarity remained unnoticed, 
mainly because of a lack of interaction between the two communities. Match- 
ing the QCD expressions with 1)1.2(1 . one discovers the hidden integrability 
properties of gauge theories [7-10]. Approximately at the same time hid- 
den integrability has been found in the M = 2 supersymmetric Yang-Mills 
theory [11]. In the three cases mentioned above, the identification goes as 
follows. 

In the so-called Bjorken kinematic limit for "hard" scattering processes, 
involving a large momentum transfer to a hadronic system, the short- 
distance perturbative QCD dynamics can often be separated (factorized) 
from the large-distance interactions and described through a set of gauge- 
invariant, local composite operators built from fundamental fields and co- 
variant derivatives. These (Wilson) operators mix under renormalization 
and their scale dependence is governed by the renormalization group (RG) 
or Callan-Symanzik equation 

H-^On{x) = Y,1nk{9)Ok{x) , (1.4) 

where ^nk is the mixing matrix given by a series in the running coupling 
constant g = g{p^)- The size of the mixing matrix is constrained by the 
symmetries and depends on the operators under consideration. The matrix 
^NK can be interpreted as a Hamiltonian acting in the space of operators 
that get mixed via the RG flow [12] with logarithm of the RG scale r = In/i 
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playing the role of the "evolution time" . In this way, the evolution equation 
(|1.4|1 takes the form of a Schrodinger equation. It turns out that for a 
certain subclass of operators, and to one-loop accuracy, the corresponding 
Hamiltonian can be identified as that of the open and/or closed Heisenberg 
magnet with the spins being the generators of the SL(2,M) group [10,13,14]. 
The number of sites in the spin chain is given by the number of fundamental 
fields involved in the composite operators and the value of the spin at each 
site is fixed by the SL{2,'R) representation to which the corresponding field 
belongs. It is different for quarks and gluons. The emergence of the SL{2, M) 
group as a symmetry group of the spin chain is not accidental, since this 
group is just a reduction of the four-dimensional conformal group 50(2,4) 
for field operators "living" on the light-cone [15,16]. Integrability allows one 
to apply the Bethe Ansatz to reconstruct the spectrum of the anomalous 
dimensions [13,14,17]. Work in this direction [18-22] has led to an almost 
complete understanding of the spectrum of anomalous dimensions of twist- 
three operators in QCD which are important for phenomenology. It has 
to be mentioned that, as a rule, QCD evolution equations only become 
integrable in the large Nc limit. There is no chance that a quantum SU{Nc) 
theory would turn out to be integrable for any Nc because the phenomena 
it describes are too complicated - nuclear forces being one of these. Indeed, 
in QCD it was found that corrections that break integrability depend on 
the quantum numbers of the operators (states) and the structure of such 
corrections are such that they lead to creation of "mass gaps" in the spectrum 
of anomalous dimensions which can be interpreted as creation of bound states 
in the corresponding quantum mechanical model [18]. 

Although historically the phenomenon of integrability was first discovered 
in QCD for operators with maximum helicity, it is, in fact, a general hidden 
symmetry of all Yang-Mills theories, which is not manifest an the classical 
level and is enhanced in its super symmetric extensions [23]. In a supersym- 
metric Yang-Mills theory, the mixing matrix in (|1.4|) gets modified due to 
the presence of additional fields. Their contribution preserves the QCD- 
type integrability and further augments it to an increasingly growing class 
of operators as one goes from pure Yang-Mills (AA = 0) to the maximally- 
supersymmetric = 4 SYM theory [24]. In particular, in the AT = 4 SYM 
theory, integrability holds in the sector of scalar operators [25,26]. This 
theory involves three complex scalars and the one-loop mixing matrix for 
the scalar operators of the type tr {<I)f^(0)<I>2^(0)<I>3^(0)} can be identified as 
the Heisenberg SO (6) spin chain with Ji + J2 + J3 sites. The SO (6) group 
is nothing but the i?— symmetry group of the model with six real scalars 
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belonging to its fundamental representation. As a natural generalization 
of these two integrable structures, discovered independently, it was found 
that the one-loop renormalization of operators involving gauge, fermions 
and scalar fields in = 4 SYM is described by a spin chain with SU (2, 2|4) 
group representations on each site [23]. 

As a second example, we consider the high-energy (Regge) asymptotics 
of scattering amplitudes. The scattering amplitude A{s, t) in QCD in the 
large- A/'c limit (here s, t are Mandelstam variables) can be written in the 
Regge limit s 3> — i as 

oo 

A{s,t)=J2-^N{s,t), AN{s,t) S^^^ , (1.5) 

N=2 

where A = g'^Nc is the 't Hooft coupling constant and the positive inte- 
ger N can be thought of as the number of (reggeized) gluons exchanged in 
the i— channel. The partial amplitudes AN{s,t) have a power-like energy 
dependence which is governed by the numbers Ef^. It turns out that the 
En coincide with the ground state energies of a completely integrable spin 
chain model [7-9]. The length of the spin chain N equals the number of the 
reggeized gluons involved, the "evolution time" is given by the total rapidity 
r = In s and the relevant degrees of freedom are the two-dimensional trans- 
verse coordinates of the reggeized gluons in the scattering plane. This model 
can be identified as a homogeneous Heisenberg magnet (|1.2() with the spin 
operators S„ being the generators of a (infinite-dimensional) unitary repre- 
sentation of the SL{2, C) group. This model can be solved by the Bethe 
Ansatz and this allows one to calculate the spectrum of Ejy [27, 28] . 

We should stress that integrability of reggeon interactions is not spe- 
cific for QCD. Since the dominant contribution to the scattering amplitudes 
comes from t— channel exchange of particles with maximal spins (gluons), it 
is not altered (in the leading logarithmic approximation) by the presence of 
additional fields in SYM theories. In other words, the phenomenon is tied 
to the gauge sector of a SYM theory. 

Finally, the third example of integrability is provided by the Seiberg- 
Witten solution for the low-energy effective action of the N = 2 SYM the- 
ory [29]. This solution is described by a Riemann surface bundled over the 
space of order parameters characterizing the vacuum state of the M = 2 
SYM theory. This surface fixes both the low-energy effective action depend- 
ing on the massless low-energy modes Uk = (tr cf)^) and the spectrum of stable 
BPS states in the theory. The Seiberg-Witten solution admits an elegant 
interpretation in the context of classical integrable models [11,30,31]. The 
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particular pattern of the integrable model depends on the matter content 

in the J\f = 2 theory, involving examples of Toda and Calogero systems, as 
well as different types of spin chains. The number of degrees of freedom in 
the underlying integrable system is equal to A^c~ 1 for gauge group SU{Nc). 
As in previous examples, the crucial point is the identification of the "time" 
variable. In the present case, the time is related to the fundamental scale of 
the J\f = 2 SYM theory, r = In Aqcd , which arises through the dimensional- 
transmutation phenomenon. The gauge-invariant order parameters turn 
out to be higher Hamiltonians of the integrable model, while the BPS spec- 
trum coincides with the action variables. 

The phenomenon of integrability has two different aspects. From the 
point of view of phenomenology, it allows one to apply the powerful Quantum 
Inverse Scattering Method [32] to solve problems which cannot be treated 
by standard techniques. Recent applications include the calculation of the 
spectrum of the multi-reggeon compound states responsible for the power- 
law rise of the scattering amplitude as a function of the energy and the exact 
solution of the evolution equations for multi-particle distribution amplitudes 
in QCD. On the theory side, integrability leaves a lot of questions to be 
answered. Currently, it has the status of an "'experimental observation" and 
its origin within the Yang-Mills dynamics remains unclear. It is natural to 
ask whether integrability survives at higher loops, what the meaning of the 
higher integrals of motion on the Yang-Mills side is and what the impact of 
non-perturbative effects is. It is unlikely that the answers can be found in 
the framework of perturbation theory. The string/gauge duality provides an 
alternative approach. It offers a unifying picture for the gauge theory both 
at weak and strong coupling regimes and allows one to relate integrability 
of the Yang-Mills theory with the hidden symmetry of a (non)critical string 
theory propagating in the curved background [33-36].^ 

Presently the gauge/string duality is rather firmly established only for the 
maximally-supersymmetric gauge theory [38]. According to this conjecture, 
at large Nc, M = 4 SYM is equivalent to a string theory with AdS^ x 5*5 target 
space and the flux of the four-form field inducing the number of colors. It is 
assumed that the radii of AdS^ and 6*5 coincide and the tension of the string 
is proportional to (g'^Nc)^/'^. At strong coupling, the tension is large and the 
string can be treated semi-classically. On the other hand, the weak coupling 
regime in the field theory is mapped into the quantum regime in the string 
sigma model. Unfortunately there exists no satisfactory quantum treatment 



Similar nonlocal Yangian type symmetries were previously discussed in the context of high-energy 
asymptotics in QCD in Ref. [37]. 
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of the string in the AdS^ x S5 background. Therefore, it is necessary to 
select proper objects both on the gauge theory and string sides to make 
their comparison explicit. 

The anomalous dimensions of the gauge theory operators appear to be 
very convenient objects for such a comparison. The dilatation operator in 
the gauge theory can be identified with the Hamiltonian on the stringy side 
so that calculation of the anomalous dimension of the operator is equivalent 
to the calculation of the energy of the corresponding string configuration 
[39,40]. Lacking the whole quantum spectrum of the string, the analysis can 
be performed only in the special, exactly solvable subsectors of the string 
theory, or for selected operators in the gauge theory which can be treated 
semi-classically on the string side. The exactly-solvable example deals with 
a string moving in the pp-background, which corresponds to the so-called 
BMN operators [39], involving scalars, in the A/" = 4 SYM theory. This is 
the only example of a gauge theory operator whose anomalous dimensions 
are known at all values of the gauge coupling constant. A more general class 
of operators with large quantum number with respect to the ii— symmetry 
group allow a semi-classical stringy treatment and the comparison of the 
stringy calculation with the perturbative expansion in Yang— Mills theory. 

The intcgrability arises along these lines in two different ways. In 
the strong coupling regime, we treat the string classically and a finite- 
dimensional integrable system of the Neumann type emerges if one assumes 
a simple ansatz for the string motion [41]. On the other hand, one can start 
with the one-loop calculation at weak coupling and consider the thcrmody- 
namical limit of the spin chain. In this limit, in the coherent state basis, the 
string sigma model becomes manifest, with each field incorporated into the 
composite operator playing the role of a single string bit [42] . 

The mapping of the integrable system to the stringy picture behind the 
low-energy eff'cctivc actions is slightly difi'crcnt. In this case, the spectral 
curve of the integrable system is part of the background since the M5 brane 
providing the world-volume for the four-dimensional theory under consider- 
ation is wrapped around it [43] . The states of the string wrapped around the 
spectral curve amount to the spectrum of BPS states in the gauge theory 
which is a counterpart of the gauge/string correspondence in this case. The 
first step towards the derivation of the stringy picture for the Regge case 
based on the semiclassical limit of the SL{2, C) spin chains was done in [44]. 

The presentation is organized as follows. In Sect. 2, we explain the phe- 
nomenon of hidden intcgrability for the scale dependence of local opera- 
tors, alias the renormalization group (RG) dilatation operator. We consider 
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three-quark operators in some detail, and using this example illustrate the 

main steps that lead to the spin-chain interpretation. This includes the 
treatment of RG equations as a Hamiltonian problem, finding an explicitly 
S'L(2)-covariant representation and identification of the conserved charges. 
The main approaches used in the quantum integrability framework are de- 
scribed, including the method of Separated Variables, the method of the 
Baxter Q— operator and semiclassical solution to the Baxter equation in 
terms of Riemann surfaces. We will demonstrate that integrability in QCD 
is restricted to operators built of fundamental fields with maximal helicity 
(and covariant derivatives), whereas the "interaction" between quarks with 
opposite helicity breaks integrability and leads to creation of mass gaps in 
the spectrum. We also give a short summary of the applications to quark- 
antiquark-gluon operators, in which case open spin chains arise. 

In Sect. 3, we consider the extension of these results to the case of 
SYM theories using the non-covariant light-cone superspace formalism due 
to Mandelstam [45] and Brink et al [46]. We give a short review of this for- 
malism and present the results of the one-loop calculation of the dilatation 
operator in the M = 2 and M = I SUSY theories. We explain how the super 
spin chains emerge in the dilatation operator, which we map to a SL(2\J\f) 
Heisenberg (super)spin chain. 

Section 4 is devoted to the high-energy asymptotics of scattering ampli- 
tudes. We explain both the differences and the similarities of this problem 
with the case of local-operator renormalization and give a short summary of 
the existing results. In this way, we identify the compound states of reggeized 
gluons in multicolor QCD as ground states of the quantum spin chain with 
SL(2,C) group. We present the exact solution to the eigenproblem for this 
integrable model based on the method of the Baxter Q— operator and discuss 
the properties of the energy spectrum within the semiclassical approach. 

Section 5 contains a discussion of the low-energy effective actions in SYM 
while in Sect. 6 we give an overview of the relation between the spectrum 
of the anomalous dimensions of composite operators in gauge theory and 
energies of the string in some background within the string/gauge duality. 
In particular, the operators with large quantum numbers will be described 
semi-classically on the stringy side. The mapping between integrable spin 
chains and strings will be briefly outlined. The flnal Sect. 7 contains a 
summary and concluding remarks. 
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2. Light-cone dominated processes in QCD 

For a practitioner, QCD as a theory of strong interactions has been mainly 
of use in scattering processes at high energies. In this case, world lines of 
quarks and gluons participating in the scattering event are close to the light 
cone. It is therefore not surprising that the separation of transverse and 
longitudinal coordinates with respect to the scattering plane proves to be 
essential. Thus, before going into details let us introduce some notation. 
Let n and n be two independent light-like vectors 

n2 = n^ = 0, n-n = l. (2.1) 

For an arbitrary four- vector, ^4^ we define 

A+ = A^n" , A_ = A^n'' , (2.2) 

and the metric tensor in the directions orthogonal to the light-cone 

= 9nv - n^n„ - nr,n^ . (2.3) 

We will also use the notation A_\_ for a generic transverse projection and 
A^^ = {0,A±,0) for a vector that only has transverse components. For 
example, 

x^" = x-v}" + x^n^ + , x^^ = g'^x„ (2.4) 

and, therefore, = 2x+X- — x\ with x\ = —x± ■ x±. 

In the parton model, hadron states are described by a bunch of partons, 
all moving in the same light-like direction, say n^, and bound in a hadron 
wave function. If, in addition to high energy, the physical process of interest 
also involves a large momentum transfer, then in a rather generic situation 
the transverse structure of the parton system is not resolved by the inter- 
action and appears to be irrelevant. For such processes, which we refer to 
as light-cone dominated, the nonperturbative hadron structure is encoded 
in suitable matrix elements of gauge-invariant nonlocal operators built of 
quark, antiquark and gluon fields located on the same light-ray (hence the 
name light-ray operators [47]) and connected by light-ray ordered gauge links 
(Wilson lines). For example, 

'il!{zin)fiPexpi^ig J du A^{un)^ tl){z2n) , 

where xp{x) is a quark field and Zj's are real numbers Zi = Xi-. In what 
follows we will not show the gauge links in order not to complicate the 



February 1, 2008 18:57 WSPC/Trim Size: 9.75in x 6. 5in for Proceedings 



BBGK 



Integrability in QCD and beyond 



275 



formulae. However, they are always implied. 

As always in a field theory, taking the asymptotic limit (here, by ap- 
proaching the light-cone) induces extra divergences in addition to conven- 
tional ultraviolet ones. Both have to be renormalized. The nonlocal light-ray 
operator can be understood as the generating function of renormalized local 
operators, e.g.. 



where D + = n- D —n- D, is the covariant derivative. Note that the 
contraction of covariant derivatives with the light-like vector picks up the 
contribution of the maximum Lorentz spin, and all arising local operators 
have the same twist = dimension - spin. The operators in 1)2. 5|) can mix with 
similar operators containing total derivatives. However, the corresponding 
mixing matrix is triangular, so that the anomalous dimensions are not af- 
fected and can be calculated by taking the matrix elements over states with 
equal momenta, in which case operators with total derivatives do not con- 
tribute. For a given A^, therefore, one is left with a single anomalous di- 
mension 7Ar, and the corresponding expression has been known since the 
pioneering work in [5]. By now it is known to three- loop accuracy. Taking 
into account the mixing with gluons makes a 2 x 2 matrix, but does not 
add complications of principle. 

Going over to operators built from three and more fields, the situation 
becomes considerably more involved. In the bulk of this section we will 
consider a particular example of the light-ray operators built out of three 
quark fields of definite helicity ^I^('L) = ^ (1 it 75) t/; at a light-like separation 



5^/ (zi, Z2, Z3) = e*^''=(/iVl)a(^in)(^^j)/3(^2n)(;iVl)7(^3n), (2.7) 



where a, /3, 7 are spinor and i, j, k color indices and the superscript refers to 
the total helicity A = 3/2 or A = 1/2. The factors /i project onto the -0+ 
components of the four-component Dirac spinor fields 



cf. Eq. (|2.5() . We will tacitly assume that the three quarks have different 





N 




^ We liope that using the same notation i/) for the quark field and the Euler i/)-function will not 
create confusion. 
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flavor. Identity of the quarks does not influence renormalization but rather 
introduces certain selection rules which pick up eigenstates with particular 
symmetries. The operators H'2.6() and H2.7() appear in the QCD theory of 
hard exclusive processes. They are used to deflne the longitudinal momen- 
tum fraction distributions of quarks inside the A-resonance and the proton, 
respectively. 

Similarly to (|2.5|) the renormalized three-quark light-ray operators are 
deflned through their Taylor expansion 

B{zi,z2,zs) = Yl '^f^f^D'^^mDi^mDi^m- (2.9) 

N ki+k2+k3=N ^' ^' 

Local three-quark operators having the same total number of derivatives all 
mix together but the ones containing total derivatives can be discarded by 
taking forward matrix elements. One is then left with a nontrivial mixing 
matrix of size -|- 1 for given A^. Eigenvalues of this matrix 'yN,q now have 
two indices: A^, which refers to the total number of derivatives, and q, which 
enumerates the anomalous dimensions (e.g., from below) for a given N. In 
the case of 53/2, we will be able to identify q with a conserved charge. 

The rest of this section is organized as follows. After a short exposition 
of the conformal symmetry in Sect. 2.1, we develop in Sect. 2.2 an approach 
to operator renormalization in which this symmetry becomes manifest. In 
this form the relation to spin-chain Hamiltonians also becomes apparent and 
this connection is elaborated in detail in Sect. 2.3. In Sect. 2.4, we use com- 
plete integrability to construct a semiclassical expansion of the anomalous 
dimensions in powers of 1/A^. In Sect. 2.5, the difference between B^^'^ 
and 5^/2 jg discussed and we find that breaking of integrability in the case 
of i?^/2 produces a mass gap in the spectrum of anomalous dimensions by a 
mechanism that can be interpreted as binding of quarks into scalar diquarks. 
Sect. 2.6 contains a short summary of other applications, in particular, the 
case of quark- antiquark-gluon operators that are related to open spin chains. 
In Sect. 2.7, we introduce the general notion of quasi-partonic operators 
which will be important for the discussion of supersymmetric extensions. 
Finally, in Sect. 2.8 we comment on the conformal symmetry breaking in 
QCD at higher orders. 

2.1. Conformal symmetry and the collinear subgroup 

Among the general coordinate transformations of the four-dimensional 
Minkowski space that conserve the interval ds^ = g^y{x)dx^^dx'^ , there are 
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transformations that change only the scale of the metric: 

a't^ui^') = ^{x)g,^u{x) (2.10) 

and, consequently, preserve angles and leave the light-cone invariant. All 
transformations belonging to this subclass form, by definition, the confor- 
mal group. It is obvious that conformal transformations correspond to a 
generalization of the usual Poincare group, since the Minkowski metric is 
not changed by translations and Lorentz rotations. Examples of specific 
conformal transformations are dilatations (global scale transformations) and 
inversion 

x^' ^ x'^" = Xxf" , x''^x'^' = ^, (2.11) 

with real A. Another important example is the so-called special conformal 
transformation 

^ ^'M = \ , , (2.12) 

1 + 2a ■ X -\- a^x^ 

which corresponds to the sequence of inversion, translation by an arbitrary 
constant vector and a further inversion with the same origin. The full 
conformal algebra in 4 dimensions includes fifteen generators: (4 transla- 
tions), M^jy (6 Lorentz rotations), D (dilatation) and (4 special confor- 
mal transformations) and is a generalization of the well-known 10-parameter 
Lie algebra of the Poincare group generated by and M^,y . 

For fields "living" on the light-ray, only a few conformal transformations 
are relevant. One of them is a special case of the conformal transformation 
in Eq. (|2.12|1 . with being a light-like vector = an^. One finds 

x_ ^ x' = — (2.13) 

" l + 2ax_ ^ ^ 

so that these transformations map a light-ray in the x_ -direction into itself. 
Together with the translations and dilatations in the same direction, x_ — > 
X- + c and x_ Ax__, they form a so-called collinear subgroup of the full 
conformal group. This subgroup is nothing but the familiar 5L(2,R) group 
of projective transformations on a line: 

2; ^ 2;' = , a5 - p-f=l, 

+ 

^z) ^ '^'{z) = (jz + 5)-2^ci> f ^ii^) , (2.14) 

\jz + d J 
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where $ is a generic quantum field (quark or gluon), a, (3, 7, S are real num- 
bers and j is the conformal spin 

j = {d + s)/2. (2.15) 

The parameter d is called the scaling dimension and it specifies the field 
transformation under the dilatations 

6d^{x) = i[^{x), n] = -{x-d + d) $(x) . (2.16) 

In a free theory (i.e., at the classical level) the scaling dimension coincides 
with the canonical dimension d'^'™, which is fixed by the requirement that 
the action of the theory is dimensionless. In the quantum theory d 7^ d^^'^, 
in general, and the difference is called the anomalous dimension. In turn, s 
is the spin projection of $ on the "+" direction: 

S_+$(z) = s$(z). (2.17) 

Here we assume that the field $ is chosen to be an eigenstate of the spin op- 
erator S |_. In the general case, one should use suitable projection operators 

to separate different spin components. 

The collinear conformal transformations are generated by the four gen- 
erators P+, M ^, D and K_ which form a subalgebra of the full confor- 
mal algebra. In order to bring the commutation relations to the standard 
form, it is convenient to introduce the following linear combinations: L_|_ = 
Li + iL2 = -«P+, L_ = Li - iL2 = (i/2)K_ and Lq = {i/2){'D + M_+). 
The action of the generators on primary quantum fields can be traded for 
the algebra of differential operators acting on the field coordinates on the 
light cone 

[L+,$(z)] = -5,$(z) = L+$(z), 

[L_, $(z)] = {z^d, + 2jz) ^z) = L_$(z) , 

[Lo, Hz)] = {zd, + j) Hz) = LoHz) , (2.18) 

where dz = d/dz. The differential operators Lj satisfy the familiar SL{2) 
commutation relations: 

[Lo, L^] = ±L^ , L+] = 2Lo . (2.19) 

The remaining generator E = (i/2)(D — M (_) counts the so-called collinear 

twist of the field $: [E, ^{z)] = ^{d - s)^{z). It commutes with all Lj and 
is not relevant for most of our discussions. 
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A local field operator ^{z) that has fixed spin projection H2.17|) on the 
light-cone is an eigenstate of the quadratic Casimir operator 

[L„ [L„ ^z)]] = j{j - iMz) = L^^iz) (2.20) 

j=0,l,2 

with the operator defined as 

L^= Yl = {Lof -Lo + L_L+, [L\Li]=0. (2.21) 

i=0,l,2 

Eqs. (|2.14l) and 1)2. 20(1 imply that the field ^{z) is transformed under the 
projective transformations according to a representation of the SL{2,M) 
group specified by the parameter j; hence the rationale of referring to it 
as the conformal spin of the field. We will see below that we are dealing 
with infinite-dimensional representations of the collinear conformal group. 
In general, different spin components of the fields have different conformal 
spin. For example, for the "good" components of the quark field, Eq. 1)2. 8() . 
s = +1/2 and the conformal spin takes the value jq = 1; for the transverse 
component of the gauge strength tensor n^Ffj^± the conformal spin equals 
jg = 3/2 and for the scalar field one has js = 1/2. 

Besides the collinear subgroup just described, one can consider another 
subgroup corresponding to transformations of the two-dimensional trans- 
verse plane = (0, xi,X2,0), introduced in ()2.4|) . This "transverse" sub- 
group is relevant for the Regge kinematics and will be considered in the 
corresponding section. 

2.2. Hamiltonian approach to operator renormalization 

In quantum theory, composite operators mix with each other under the 
renormalization group flow in the cut-off parameter. This flow is driven 
by the dilatation operator, 

[0,D] = i (^d--^ - /^l; - /5(5)^) O , (2.22) 

where dean is the canonical dimension of the composite operator O. This 
equality is understood as an operator insertion into a Green function and is 
known as the Callan-Symanzik equation. Conformal symmetry of the QCD 
Lagrangian (for massless quarks) leads to important constraints on the above 
operator mixing. Indeed, the one-loop renormalization is governed by the 
infinite parts of the one-loop Feynman diagrams which have the symmetry 
of the classical Lagrangian. It follows that the QCD operators belonging to 
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^--.-^----^ ^----^----^ ^----^----^ 

a b c 



Figure 1. Examples of a "vertex" correction (a), "exchange" diagram (b) and self-energy 
insertion (c) contributing to the renormalization of three-quark operators in Feynman 
gauge. Path-ordered gauge factors are shown by the dashed lines. The set of all diagrams 
includes possible permutations. 

different representations of the (collinear) conformal group cannot mix under 
renormalization at leading-logarithmic accuracy. This simplification proves 
to be crucial, e.g., for the QCD description of hadron form factors at very 
large momentum transfers. The aim of this section is, first, to explain this 
symmetry of the QCD renormalization group equations using a simple ex- 
ample and, second, to re-interpret these equations as a quantum-mechanical 
Hamiltonian problem. 

To make the conformal symmetry explicit, it is convenient to work with 
light-ray operators 1)2. 6|) . (|2.7() directly, rather than decompose them in local 
operators. The renormalization group equations for light-ray operators can 
be written as [47-49] 

!^f^^ + P{g)^^B = M-B, (2.23) 

where H is a certain integral operator corresponding, to the one-loop accu- 
racy, to contributions of the Feynman diagrams shown in Fig. 

To simplify notations we factor out the QCD coupling constant, the color 
factors and trivial contributions of the self-energy insertions: 

(2.24) 

where Cp = {N^ — l)/{2Nc) is the Casimir operator in the fundamental rep- 



The color factors of the diagrams shown in Figs. and b can be calculated using the identity 
^i' i' k^l' i^'ji j — ~^ijk{^ + l/^c)/2 with t" being the generators of the SU{Nc) in the quark 
representation. 



9 

8^2 



(1 + 1/Nc)n + 3Cf/2 
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resentation of SU{Nc). In the Feynman gauge, the gluon exchange diagram 
in Fig. vanishes unless the participating quarks have opposite hehcity. 
The renormahzation of the A = 3/2 operator -B3/2 l|2.(ij) is therefore deter- 
mined by the vertex correction in Fig. ^ alone. By explicit calculation one 
finds [50, 51] 

W3/2 = + W23 + ^13 > (2.25) 

where TC^/^ are the two-particle kernels involving the i-th and k-th quarks, 
for example, 

Jo a L 

+ a[Bizi,z^^,Z3)-B{zi,Z2,Z3)]], (2.26) 

with a = 1 — a and z^j^ = ZiOt + z^a. 

In the case of -B1/2, the vertex correction remains the same, but one has 
to add the contributions of gluon exchange between the quarks with opposite 
helicity. One obtains 

7^1/2 = W3/2 - ^12 - ^23 1 (2.27) 

where we assume that the first and the third quark have the same helicity, as 
in Eq. 1)2. 7|) . The kernels Hff^ act on i-th and A;-th arguments of the nonlocal 
operators only, and can be written in the form 

B{z,) = I VaB{z1^,z^l,Z3) , (2.28) 

with the integration measure Da defined as 

J'Da = J dai (iQ2 das 5(1 — «! — 02 — 03) • (2.29) 

The expected conformal invariance of the evolution equation for bary- 
onic operators implies that the two-particle kernels Tiik commute with the 
generators of the SL{2) transformations La defined in (|2.18|) . To show this, 
consider the following expression that generalizes both 1)2. 26() and (|2.28l) : 

ni2B{zi,Z2,Z3) = j 'Dauj{ai,a2)B{zi - aiZi2,Z2 + a2Zi2,Z3), (2.30) 

where Z12 = zi — Z2- The integration measure is defined in 1)2. 29() . This 
operator has a simple meaning: acting on the three-particle nonlocal oper- 
ator B{zi, Z2, Z3) it displaces the quarks with the coordinates zi and Z2 on 
the light-cone in the direction of each other. The vertex correction in 1)2. 26() 
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is obtained with uj'"{ai,a2) = ^ (^§7§f j the exchange contribution in 
(|2.26|) corresponds to (01,02) = 1. 

It is easy to see that, for this ansatz, [Hi2,L^] = ["^12,-^^0] = for an 
arbitrary function uj(ai,a2), whereas the condition [Hi2,L-] = leads to 
the following constraint: 

^^oioi + 2aiji^ uj{ai,a2) = ^^^020:2 + 202^2^ ^^(01,02) , (2.31) 

where a = 1 — a. Its general solution has the form 

/ ^ -2ii-2-2i2-2 f ai<^2\ fn on\ 

uj{ai,a2) = a,-"- a, 9^ _ _ , (2.32) 

\a1a2J 

with an arbitrary ip. Note that the conformal spins jk = 1 for all the 
three quark fields entering (|2.6j) and (|2.7|) so that the prefactor q;^''^~^q;2''^~^ 
in (|2.32|) disappears and the conformal symmetry requires that the kernel 
io{ai,a2) is an arbitrary function of the ratio which is indeed the case 
for both expressions ((T^ and 

Once conformal symmetry of the two-particle kernels is established, the 
group theory tells us that Tiik may only depend on the corresponding two- 
particle Casimir operators L^^. To find the functional form of this depen- 
dence one has to compare their action on a suitable basis of trial functions. 
For definiteness, let us find 7ii2 as a function of L^2- To this end, it is enough 
to compare their action on the homogeneous polynomials of two variables zi 
and Z2- 

B{zi,Z2,Z3) > bn{zi,Z2) 

which we choose to be eigenfunctions of the operator 

= -did2zl2; zi2 = zi- Z2. (2.33) 

It is easy to see that the polynomials defined thus form an infinite- 
dimensional representation of the SL{2) group on which the operators 
L-|_ = + L2,+ and L_ = -Li,- + L2- act as the lowering and raising 
operators respectively. It is thus sufficient to consider only the functions 
(polynomials) annihilated by L_|_, or equivalently, the highest weight of the 
representation, since all other eigenfunctions of can then be obtained by 
a repeated application of L_. Since = —{di + 82), the latter condition 
is simply translation invariance, which leaves one with 

6(zi, 22) = (21-^2)" = ^1^2 , n = 0,1,2,... (2.34) 
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An explicit calculation gives 

^?2^r2 = (™ + 2)(n + lK2> 
WM = 2[V'(n + 2)-V(2)]zr2, 

nt,z^, = l/[{n + 2){n + l)]z^„ (2.35) 

where ip{x) = d\nT{x)/dx is the Euler function. To cast (|2.35() in an 
operator form, define J12 as a formal solution of the operator relation 

LI2 = Jui'Ju - 1). (2.36) 

The eigenvalues of J12 equal ji2 = n + 2 and specify the possible values of 
the sum of two j = 1 conformal spins of quarks in the (12)-pair, cf. (|2.21j) . 
Then 

Wi2 = 2[V(Ji2)-V'(2)], 

W12 = l/[Ji2(Ji2 - 1)] = 1/L?2 • (2.37) 

Substituting the representation 1)2. 37() into ()2.25|) and ()2.27p one obtains a 
Schrodinger equation for three particles with the coordinates zi, Z2 and 23 on 
the light-ray line. The 'Hamiltonians' 'H3/2 and ^1/2 entering this equation 
for different baryon states have a pairwise structure and are expressed in 
terms of the corresponding two-particle Casimir operators 1)2. 37() . 

The precise form of the Schrodinger equation and the meaning of the 
corresponding "wave function" depend on the representation that is used 
for the SL{2) generators. Looking for polynomial solutions in light-cone 
coordinates $(zi, 22,-23) as in ()2.34|) is one possibility: 

n ■ ^N,q = £N,q'^N,q ■ (2.38) 

Here N refers to the degree of the polynomial alias the total number of 
derivatives, and q enumerates the energy levels. 

Another option is to go over to local operators. A generic local operator 
with derivatives can be written as the sum of monomials entering the 
expansion (|2.9() with arbitrary coefficients 

o= Yl ckM^Di^mDi'mD^'m, (2.39) 

ki+k2+k3=N 

and can be represented by a polynomial in three variables 

^{xi,X2,X3)= Ck.k^k^Xi'x^^x'^' . (2.40) 

ki+k2+k3=N 
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Note that ^{xi) serves as a projector, separating out the contribution of 
the local operator Oiji to the nonlocal operator B{zi), which can be made 
explicit by writing 

= ^{di,d2, d-i)B{z,,Z2, Z3)U=0- (2.41) 

Therefore, one can call ^{xi) the coefficient function of a local operator. 
Since a local operator is completely determined by its coefficient function, 
diagonalization of the mixing matrix for operators can be reformulated as 
diagonalization of the mixing matrix for the coefficient functions. Requiring 
that Oq/ 1)2.41(1 is multiplicatively renormalized, one ends up with a matrix 
equation in the space of homogeneous polynomials of degree of three 
variables 

n ■ ^N,g = £N,q^N,q , (2.42) 

where the Hamiltonian 7i is given by the same expression, Eq. (|2.37|) . in 
terms of the two-particle Casimir operators, but the SL{2) generators have 
to be taken in a different, adjoint representation (see [18,22] for details). The 
eigenvalues of Ti in (|2.38|) and l|2.42|l are of course the same; they correspond 
to the anomalous dimensions 



a. 



^^'^ = 2^ 



[l + l/N,)£N,g + S/2CF ■ (2.43) 



The wave functions in (|2.38() and (|2.42|) are related by a certain integral 
transformation, see [18]. The advantage of using ^(zi) is that is this rep- 
resentation the highest weight vector condition L+$(zj) = reduces to the 
requirement of the translation invariance of ^{zi) = 0. On the contrary, 
the corresponding equation for the operator coefficient functions is more 
complicated (see also [52]). 

The Hamiltonian in (|2.42() is hermitian with respect to the conformal 
scalar product [16]. This implies that the anomalous dimensions take real 
quantized values and the corresponding eigenfunctions are mutually orthog- 
onal with the weight function xi 3:2X3 

j I)XXIX2X3 'ifN,q{Xi)^N,q'{Xi) ~ 6g^g' . (2.44) 



The construction that we have described was first suggested in [12]. It is 
general and can be used in all other cases. 
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2.3. Complete integrability and non-compact Heisenberg 
magnets 

The Hamiltonian ^3/2 possesses an additional "hidden" symmetry. One can 
construct an integral of motion (conserved charge) that commutes with ^3/2 
and with the SL[2) generators: 

[Q,L,] = [2,^3/2] =0, (2.45) 

where the explicit expression in the first line refers to the representation of 
the generators in Eq. (|2.18|) . The existence of a non-trivial integral of motion 
makes the corresponding Schrodinger equation (|2.38j) completely integrable. 
Indeed, Eq. (|2.45j) tells us that W3/2, Q, the total three-quark conformal spin 
= Z]i=o,i,2(-^i,i + ^2,i + Ls^if and its projection Lq = Li,o + ^2,0 + ^3,0 
can be diagonalized simultaneously. Since there are only three independent 
degrees of freedom, this implies that the Hamiltonian "^3/2 is a (compli- 
cated) function of the conserved charges. As a reward, one can replace the 
Schrodinger equation by a simpler problem of diagonalization of the con- 
served charge Q. The eigenvalues of Q provide one with a "hidden" quantum 
number that can be used to parameterize the spectrum of the Hamiltonian, 

S3/2=£3/2iN,q).'' 

It turns out that ^3/2 can be viewed as the Hamiltonian of the non- 
compact SL{2,M) Heisenberg spin magnet. As a hint, notice that the ex- 
pression for the Hamiltonian "^3/2 as the sum of the Euler functions, 
Eqs. ()2.25|) and 1)2. 37() . resembles the definition of generalized SU{2) Heisen- 
berg magnet, Eqs. ()1.2|) and (|1.3|) . The important difference between the 
two expressions is, however, that the two-particle "spin" J12 in (|2.37|) takes 
an infinite set of quantized values J12 = n + 2 with n > 0, while for the same 
operator in ()1.2() one has < J12 < 2s. In other words, in these two cases 
we deal with the mfinite-dimensional representations of the S'L(2,R) group 
and iinite-dimensional representation of the SU{2) group, respectively. 

To establish the exact correspondence, we can think of three quark op- 
erators inside B{zi, Z2, Z3) as three sites of the lattice and reinterpret the 
generators of the collinear SL{2,M) subgroup, Lo,L+ and L_, acting on 



For the local conformal operator with the coefficient function \E'(xi) that satisfies the highest 
weight condition L+<E'(a;i) = one finds Lo*{a;i) = J^{xi) and L'^^{xi) = J{J - l)'<!>{xi). Here 
J = N + 3j is the total conformal spin with j = I being the conformal spin of the quark field. 
Hence the dependence of the energy eigenvalues on the conformal spin (and spin projection) can 
be traded for the dependence on A'^, the total number of derivatives. 



1,2008 18:57 WSPC/Trim Size: 9.75in x 6. 5in for Proceedings 



BBGK 



286 A. Belitsky, V. Braun, A. Gorsky and G. Korchemsky 

A;— th quark field as defining the spin operators in k—th site 

S+ = zldk-2szk, S; = -dk, Sl = zkdk-s, (2.46) 

with = i: iS"^ and spin s related to the conformal spin of the quark 
field s = —j = —1. 

The spin operators (|2.46|) are the generators of a unitary representation of 
the SL{2, M) group of the discrete series. The definition of the integrable spin 
chain is based on the existence of a fundamental operator, the i?— matrix, 
Rkm{u), which acts on the coordinates in sites k and m, depends on an 
arbitrary complex parameter u and satisfies the Yang-Baxter equation [32] 

Ri2{u)Ri^{v)R2z{u -v) = R2z{u - v)Rn{v)Ri2{u) , (2.47) 

with u and v being arbitrary complex spectral parameters. For infinite- 
dimensional representations of the SL{2, M) group of the discrete series, the 
solution to Eq. 1)2. 47() is given by [3] 

RUu) = i-iV'^- l^'^;"' l':"] , (2.48) 

r( Jfcm + tu) 

where J^m is the two-particle spin operator, cf. Eq. (|2.36|) . The Hamiltonian 
of the SL{2,M) Heisenberg magnet of length L is defined in terms of the 
i?— matrix as [4] 



n=l 



(2.49) 

u=0 



with Hl^l^i = Hl i. It is easy to verify that for L = 3 the Hamiltonian 



(|2.49j) coincides with "^3/2 up to an irrelevant additive constant. 

Complete integr ability of the Hamiltonian "^3/2 allows one to solve the 
Schrodinger equation by applying the Bethe Ansatz. The exact expression 
for the energy spectrum is given by 



''^ du Q{u — ij) 



(2.50) 

u=0 



where the function Q{u) is the eigenvalue of the Baxter operator for the 
SL{2,M) magnet [53]. It satisfies the finite-difference Baxter equation 

t{u)Q{u) = {u + ijfQ{u + i) + {u- ijfQ{u - i) (2.51) 

with the transfer matrix 

t{u) = 2u^ -{N + 3j)iN + 3j -l)u + q (2.52) 
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and j = 1 being the conformal spin of the quark field, under the addi- 
tional condition that Q{u) is polynomial in the spectral parameter. These 
conditions determine Q{u) up to an overall normalization and allow one to 
establish the quantization conditions for the charge q. The polynomial Q{u) 
can be parameterized by its roots 



with non-negative N related to the total SL{2, M) spin of the magnet J = 
N + 3j. Substituting this expression into the Baxter equation H2.51() . one 
finds that the roots Xj satisfy the Bethe equations for spin s = —1. The fact 
that the spin is negative leads to a number of drastic differences compared 
to conventional "compact" SU{2) magnets. In particular, the Bethe roots 
take real values only and the number of solutions is infinite [9] . 

Solving the Baxter equation H2.51|l for = 0, 1, . . . one can reconstruct 
the spectrum of the Hamiltonian "^3/2 and, as a consequence, determine the 
exact spectrum of the anomalous dimensions of the maximal helicity baryon 
operators, see Fig. |21 

We would like to stress that the function Q{u) has profound physical 
meaning: It defines the eigenstates of the <S'L(2,M) magnet in the represen- 
tation of the Separated Variables [54] . Due to complete integrability of the 
Hamiltonian ^3/2, its eigenfunctions '^N,qizi, Z2, Z3) have to diagonalize the 
conserved charge Q defined in (|2.45|) . In general, the solutions to the differen- 
tial equation Q'$N,q{zi, Z2, Z3) = g^'iv,g(^i, ^^2, ^3) are complicated functions 
of the light-cone variables Zk- Obviously, the spectrum of "^3/2 and Q does 
not depend on the representation in which these operators are defined. The 
main idea behind the method of Separated Variables is to perform a unitary 
transformation on the Hilbert space of "^3/2 



W3/2^Wf/f ) = Uns/2U^ , ^N,g ^ = U^N,g (2.54) 



in order to go over to a representation in which the eigenfunction ^j^^ 
is factorized into a product of functions depending on a single variable. 
Denoting the coordinates in the SoV representation as 



N 



Q{u) = ll{u-Xj), 



(2.53) 



(^^1, 22, 2:3) 



SoV 



{zo,Xi,X2) 



(2.55) 



with zq being the center-of-mass coordinate of the system of three quarks on 
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the light-cone and xi^2 being new cohective coordinates, one has 

<°/^=e*^^''Q(xi)g(x2). (2.56) 

The exphcit form of the unitary transformation to the Separated Variables, 
Eqs. ()2.54l) - H2.56() . was found in [55]. Here p is the total (light-cone) 
momentum of three quarks Ylk=i''''^Zk^N,q{zi, Z2, Z3) = N,qizi, Z2, Z3), or 
equivalently the "— " component of the total momentum of the SL{2,M) spin 
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chain. The function Q{x) defines the wave function in the separated vari- 
ables, Eq. (|2.56|) . In this way the Schrodinger equation for 'ifN,q is translated 
into the Baxter relation for Q{x), Eq. (|2.51j) . We conclude that the Baxter 
equation encodes the complete information about the system - its solutions 
Q{x) govern the energy spectrum, Eq. (|2.5U|) . and define the wave functions 
in the representation of Separated Variables, Eq. 1)2.56(1 . 

2.4. Semiclassical expansion 

The spectrum of the anomalous dimensions shown in Fig. |21 exhibits remark- 
able regularity. One can reveal its origin by making use of the interpre- 
tation of Q{u) as the wave function in the separated variable. Assuming 
a parametrization Q{x) = exp(^S(a;)), one can solve the Baxter equation 
1(2. 51(1 by making a semiclassical expansion for h ^ 0. Along this line of 
reasoning, the role of the Planck constant is played by the (inverse) total 
conformal spin 1/J = 1/{N + 3j) for ^ 1. 

To be more precise, we introduce a small parameter rj 

77 = [J(J-l)]-i/2^ J = N + 3j (2.57) 

and replace u = x/r] as the parameter in Eq. ((2. 51(1 . Introducing the notation 
f{x) for the rescaled Q-function 

fix) = i^\x/r,fQ{x/r,) , q = qrf'. (2.58) 

the Baxter equation ((2.511) takes the form 

fix + ir]) + f{x - irj) = x~^t(x)/(x) , (2.59) 

where in the expression ((2.52(1 for the transfer matrix t{x), the parameter 
q is replaced by q defined in ((2.58() . We now look for the solutions of the 
rescaled Baxter equation in the form of the WKB expansion 

fix) = exp (^^Soix) + iSiix) + , (2.60) 

where each term is assumed to be uniformly bounded and the expansion is 
assumed to be convergent. Substituting the ansatz ((2.60(1 in Eq. ((2.59(1 and 
expanding both sides to order -q, one obtains [17] 

2coshS'^(x) = tix)/x^ , S[ix) = ^Sl^ix)cotliS'oix) , (2.61) 
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where a prime denotes a derivative with respect to x. To this accuracy, the 
solution of the Baxter equation (|2.59|) is given by 

As usual, the leading term 5*0(2;) defines the classical action function of the 
underlying quantum system. Introducing the notation 

Li!{x) = exp(5o(x)) = exp{p{x)) (2.63) 

one can rewrite the first equation in (|2.61l) as 

1 -3 

UJ + - = x^h{x) . (2.64) 

This equation establishes the relation between the classical momentum p and 
the coordinate x for "equal energy levels" defined by the charge q. Taking 
into account that uj is real on classical trajectories, we obtain the conditions 
that define the classically allowed regions of motion for quarks as 

t\xj)>4x^j, (j = l,2), (2.65) 

where the subscript refers to the j— th allowed band. 

Next, we continue the relation (|2.64j) into the complex domain, define 
y = x^{ol> - ^) and rewrite fU^^ as [17] 

T : =t^{x) -4x^ = {-x + q){Ax^ -x + q), (2.66) 

where t{x) = 2x^ — x + q is a polynomial depending on the charge q. The 
curve ()2.66() determines an elliptic Riemann surface T equipped with a mero- 
morphic differential cISq = p{x)dx defined in 1)2. 61() . Integrating this differen- 
tial along some path that terminates at x one obtains the "action" function 
dSo = p{x)dx. In general, for arbitrary values of the conserved charge 
q, this function is single- valued on the Riemann surface but double- valued 
on the complex x— plane. We recall that the same "action" function defines 
the wave function in the separated coordinates, Eq. 1)2.62(1 . which has to be 
a single-valued function on the real x— axis. To satisfy this requirement, the 
charge q has to satisfy certain quantization conditions. In the semiclassical 
approach, they take the form of standard Bohr-Sommerfeld quantization 
conditions imposed on the solution in Eq. ()2.62j) 

ak = dSo = 2^ri {rik + ^ + 0{rf) (2.67) 
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with k = 1,2 and being integer, such that ni + n2 = N. Here the inte- 
gration goes over the closed contours on the Riemann surface encirchng 
the intervals of the classical motion on the real x— axis, Eq. H2.65() . Solving 
H2.67(l . one finds the spectrum of quantized charges 



with 7] defined in (|2.57|) . The integer ni parameterizes different "trajectories" 
as shown in Fig. [2 Note that all non-zero eigenvalues of Q come in pairs: 
If q is an eigenvalue, then —q is also an eigenvalue. As a consequence of 
the symmetry of quantization conditions in the separated coordinates, the 
trajectories q{N, rii) have the reflection symmetry q{N, ni) = —q{N, N—ni), 
which maps positive values of q on the ni-th trajectory into the negative q 
on the 77.2 = (A^ — ni)-th trajectory. 

Finally, the energy eigenvalues £]\f,q can be calculated as a function of the 
spin N and the charge q with the help of Eq. 1)2. 5U() . At the semiclassical 
level, one obtains 



%q = 2 In 2 - 6 + 67E + 2 3f?e J] V(l + i^k/v^) + 0{ri^) (2.69) 



where 6k are the roots of the transfer matrix t3{6k) = 26^ — Sk + q = 0. Com- 
bined together, Eqs. (|2.68|) and (|2.69|) provide a good quantitative descrip- 
tion of the spectrum of the anomalous dimensions for helicity— 3/2 baryon 
operators. Note that £N,q = ^N-q so that the anomalous dimensions are 
double degenerate, apart from the lowest eigenvalues for even N which cor- 
respond to q = 0. The origin of this degeneracy can be traced to the permu- 
tation symmetry between the three quarks [18]. 

2.5. Breakdown of integrability 

There is no reason to expect that QCD is an integrable theory, simply be- 
cause of the complexity of phenomena that it is supposed to describe. This 
complexity is certainly related to breaking of the conformal symmetry, but 
even so breaking of integrability in cases where the conformal symmetry is 
preserved may lead to important clues. The renormalization of the three- 
quark operator 5^/^ provides one with a very instructive example of what 
can happen. As mentioned above, the operators B^/"^ (integrable) and B^/"^ 
(non-integrable) are related to quark distributions at small transverse sep- 
arations inside the A-isobar and the nucleon, respectively, and it is natural 
to speculate that the difference in light-cone dynamics in these two hadrons 



q{N,ni) 



\/27 . 




-g) +...]. (2.68) 



3 



k=l 
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Figure 3. The flow of energy eigenvalues for the Hamiltonian 'H(e) for N = 30 (see text). 
The solid and the dash-dotted curves show the parity-even and parity-odd levels, respec- 
tively. The two vertical dashed lines indicate 7^3/2 = 7i(e = 0) and 7^i/2 = Ti-ie = 1), 
respectively. The horizontal dotted line shows position of the "ground state" that corre- 
sponds to g = 0. 

is related to the loss of integrability in the latter case. Let us study this 
difference in some detail. 

The scale dependence of B^^'^ is driven by the Hamiltonian TC1/2 defined 
in 1)2. 27() . which differs from TCs/2 by the two terms corresponding to gluon 
exchange between quarks of opposite chirality (see Fig. P): 

= W3/2 + V, V = -(-^ + ^]. (2.70) 

V^12 ^23/ 

To visualize the effects of this additional interaction it proves to be conve- 
nient to consider a somewhat more general Hamiltonian 

n{e)=ny2 + eV, (2.71) 

with e being a new coupling constant. Ti.{€ = 0) reproduces "^3/2 whereas 
Ti.{e = 1) coincides with the Hamiltonian Ti.i/2- Thus, the spectra of 7^3/2 and 
Tli/2 are related to each other through the flow of the energy levels of Ti.{e) 
from e = to e = 1 (see Fig. Note that for e 7^ 0, the cyclic permutation 
symmetry between the three quarks is also lost, but 'H(e) is still invariant 
under permutations of the first and the third quarks; ['H{e),Vi3] = 0. As 
a result, the degeneracy between parity-odd and parity-even eigenstates is 
lifted and, in fact, the flows of levels with different parity are completely 
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independent from one another. 

The spectra in Fig. exhibit the following characteristic features: 

- In the upper part of the spectrum the effect of e— proportional terms 
on the spectrum of the 'unperturbed' Hamiltonian TC{€ = 0) = "^3/2 
is very mild. While at e = the energy levels are double degenerate, 
their splitting at e 7^ remains (exponentially, as one can show [18]) 
small for large N. 

- For e > 0, the two lowest levels are decoupled from the rest of the 
spectrum and fall off with e almost linearly. In the large— limit 
they become separated from the rest of the spectrum by a finite "mass 



This structure suggests that the difference eV = 7i{e) — Ti.{e = 0) can be 
considered as a perturbation for most of the levels, but not for the few lowest 
ones (for large N). Moreover, for e > 0, which corresponds to attractive 
interaction, we seemingly witness creation of a certain bound state. 

To formalize the argument, one has to evaluate the matrix elements of 
V between the eigenstates of "^3/2 find compare them with the energy split- 
tings for the "unperturbed" Hamiltonian. By explicit calculation one con- 
firms that the perturbation theory in V is justified for large for the upper 
part of the spectrum, while several lowest energy eigenstates are affected 
strongly and have to be re-diagonalized. It turns out that the matrix ele- 
ments for low-lying energy eigenstates can be computed in terms 
of eigenvalues of the cyclic permutation operator 



'P^'g(xi,X2,X3) = ^'g(x3,Xi,X2) = 6 *g (xi , X2 , X3) 



(2.72) 



where the phases (j)q take quantized values 4>q = 0, ^ (since = 1). 
One obtains [18] 



TT' 



.2 



1 



9 IniV 



COs(0g - (j)q>) 



(2.73) 



It is easy to see that the possible values of cos((/)g ± (pgi) are 1 and —1/2 
depending on whether the phases (pq and (f>qi coincide. Introducing an integer 
k = 0, ±1, ±2, ... to enumerate quantized values of q starting from the ones 
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with the lowest absolute value, we get 



{^g'\eV\^q) = -gAk'k, ^k'k 



( 1 



_ 1 

2 

2-^22 

1 _i 1 _1 

2 2^2 

1 _1 

2 2 



1 



1 



1 



V 



(2.74) 



where the dependence on e and r] is absorbed in the 'effective coupling' 
Notice the ~ 1/ In falloff for large N , which is slower than the 



9 



9in TV- 



level splitting for "^3/2! ^-£■3/2 ~ ) confirming our expectation that 

~ InA^ lowest levels have to be re-diagonalized. Neglecting, for a moment, 
the splitting between unperturbed levels up to some number /cmax ~ InA^, 
the true eigenstates would coincide with those of A^k' (of the size /cmax); 
and the corresponding eigenvalues would define the energies. Remarkably, 
the spectrum of Kky is extremely simple. One can easily convince oneself 
that h-kk' has two and only two nonzero eigenvalues which are both equal to 
fcmax- Remembering that g ^ \/\d.N and Ar^ax ~ In A^, this implies that two 
energy levels will get shifted by the finite amount (^^max ~ 

0{\n^N) while 

all the other ones remain exactly degenerate to this accuracy. Since (|2.74j) 
was derived up to corrections of order ~ 1/ln^ N , this implies that the true 
energy shift for all levels apart from the lowest two is, at most, ~ 1/ln^ N . 
This explains the creation of the mass gap observed in Fig. |31 

The size of the gap can be calculated for small e by constructing an 
effective Hamiltonian for the low-lying eigenstates [18]. The idea of the 
solution is to interpret the integer /c as a discrete momentum variable. The 
corresponding wave functions in configuration space correspond to Bloch- 
Floquet waves and the resulting Schrodinger equation turns out to be a 
generalization of the famous Kroning-Penney model of a single particle in a 
periodic (5-function potential. The result for the mass gap reads 



A(e) = Abound — Eq 



72C(3) 



(2.75) 



Extension of this approach to the physically interesting case e ~ 1 is non- 
trivial. However, the mass gap can easily be calculated in the opposite limit 
e ^ 1, and the result for e = 1 found by matching the two expansions using 
Fade approximants: A1/2 = A(e = 1) = -(0.32 ± 0.02) [18]. 

Creation of the mass gap implies that the wave function of the ground 
state is drastically modified. For the integrable case, the eigenfunction cor- 
responding to the lowest energy and hence q = can be found exactly. The 
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physical interpretation of this result is most transparent in coordinate space. 
Up to contribution of operators with total derivatives, one can represent the 
'ground state' of three quarks with the same helicity in a concise form as the 
nonlocal light-cone operator 



The Taylor expansion of (|2.76|) at short distances, 212,-232 0, generates 
the series of local multiplicatively renormalizable three-quark operators with 
the lowest anomalous dimension for each even A^. Note the integration in 
Eq. H2.76() with unit weight over the position of the quark in the middle 
that goes in between the light-cone positions of the other two quarks, up to 
permutations. If renormalization of the operator is interpreted as interac- 
tion, this unit weight can in turn be interpreted as the statement that the 
quark in the middle is effectively 'free': In the 'ground state' with the lowest 
"energy" , the interaction of the quark in the middle with its right and left 
neighbors exactly compensate each other. 

For the B^^"^, case the ground wave function changes dramatically and in 
the limit of large becomes, in the same sense as ()2.76|) . 



where the two solutions correspond to states with different parity (cf. Fig.|21). 
Formation of the mass gap in the spectrum of anomalous dimensions is, 
therefore, naturally interpreted as due to binding of the quarks with opposite 
helicity into scalar diquarks. 

Note that while the expression ()2.76|) for the eigenfunction is exact, the re- 
sult 1)2. 77(1 is only valid in the asymptotic In — > 00 limit. In the coordinate 
space picture, the restriction to large N is translated to the condition that 
the light-cone separation between the same helicity quarks is very large to 
allow for the formation of a diquark. In other words, the diquark itself has a 
large size. In the momentum space, the result means that at sufficiently large 
normalization scale jj? the quark carrying a very large momentum fraction 
is more often with the same helicity as of the parent baryon. This observa- 
tion seems to be in qualitative agreement with phenomenological models of 
baryon distribution amplitudes. 



b: 



,(9=0) 

'3/2 




a, 6=1, 2, 3 



(2.76) 



± M)(«3)5(a2 - 03)1 , (2.77) 
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2.6. Open spin chains 

There are other multiparticle operators which represent phenomenological 
interest and which were studied along the same hne as the three-quark opera- 
tors. For the twist-three sector, in addition to three-quark operators that we 
have addressed above one considers quark-antiquark-gluon and three-gluon 
operators. They arise in particular in connection with the structure function 
g2{x,Q'^) measured in polarized deep- inelastic scattering. This is attracting 
increasing interest as it provides a direct measurement of three-parton corre- 
lations in the proton. The scale-dependence of §2(30, Q^) can be traced to the 
renormalization-group equations for the quark-antiquark-gluon operator [56] 

S^{zi,Z2,Z3) = ^^{zi)[igh+{z2) ± 9F±+{z2h5]Mz3) , (2.78) 

where F^i, = \e^yp\F°-'P^t°- , F^,y = Fj^^t"'. In the flavor singlet sector, one 
also has to include the three-gluon operator 

d{zuz2,z^) = '-gr''FUzi)FlA^2)FU^s) , (2.79) 

which mixes with the C-even part of the flavor singlet quark-antiquark-gluon 
operator operator (|2.78j) . 

To begin with, consider the flavor non-singlet sector. In this case the 
three-gluon operator drops out and the evolution of simplifies dramat- 
ically in the limit of large number of colors, Nc 00. The corresponding 
Hamiltonian reads, e.g., for [10,12] 



a.. 



* 27r 



VgoiJu) + UGg{J23) + O {l/Nc) , (2.80) 



where 



v,g{j) = ^ ( J + i) + V ( J - i) - 2V(i) - 1 , 

UcgiJ) = V (^ + ^) + V' (J - i) - 2V(1) - I , (2.81) 

and it turns out to be equivalent to the Hamiltonian of the completely inte- 
grable open Heisenberg magnet [10]. The corresponding conserved charge is 
equal to 

Qs+ = {^125-^23} ~ 2^12 - 2^23 ' (2.82) 

where {,} stands for an anticommutator. Properties of the Hamiltonian 
H2.80() were studied in detail in Refs. [10, 13, 14, 19] and the corrections 
0{1/N^) to the spectrum were calculated in Ref. [21]. The most interesting 
result is that the lowest energy, albeit the lowest anomalous dimension in 
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the spectrum of quark-gluon operators can be exactly found, and the cor- 
responding eigenfunction coincides to leading logarithmic accuracy with the 
structure function §2 itself. Thus, a simple evolution equation arises [57]. 

The study of the flavor singlet sector is much more challenging [22] . Ignor- 
ing quarks for a moment, the evolution kernel for the three-gluon operators 
1)2. 79() can be mapped to the Hamiltonian of closed Heisenberg magnet with 
three "gluonic" sites of the spin jg = 3/2, whose integrability is "softly" 
broken by additional terms [13].^ This case was studied using the same 
techniques as the three-quark operators, although the algebra becomes more 
complicated because of the higher spins involved. The most interesting part 
appears to be the mixing between quark-gluon and three-gluon operators 
which for large can be interpreted as describing the interaction between 
open and closed Heisenberg magnets. It turns out that this mixing has 
rather peculiar features, which we cannot discuss in detail in this review, 
but the outcome is that the mixing can to a large extent be reduced to a few 
participating levels. Identifying important degrees of freedom (at least for 
large N) one can write down an approximate two-channel evolution equation 
for the structure function g2 in terms of transverse spin densities [22] . 

2.7. Quasipartonic and non-quasipartonic operators 

The operators we have discussed above all fall into the same class of the 
so-called quasipartonic operators, introduced in Ref. [12]. Since an offspring 
of this formalism will be useful in the discussion of supersymmetric exten- 
sions of QCD, let us elaborate upon it. The approach heavily relies on the 
light-cone gauge, = 0. In this gauge, only physical degrees of freedom 
propagate, i.e., "good" components of the quark field V'+ and circularly po- 
larized gluons . Thus the identification of the elementary fields entering 
the Lagrangian with partonic states is straightforward - hence the name 
quasipartonic operators. Introducing the helicity operator as^ 

h = e^e^^ , ei = (ei)* = ^(1,^) , (2.83) 

one decomposes all fields in components with a fixed helicity. In particular, 
the left V+ and right quarks, together with (anti)holomorphic compo- 
nents of the two-component propagating gauge potential A = -i=(A^ + iA'^) 



^ Tlie exactly integrable chain arises for maximal-helicity gluonic operators [13] . 
' Do not confuse helicity h with the light-cone spin, discussed above. 
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and A = A* are eigenfunctions of the helicity operator corresponding to 

hA{x) = A{x) , hA{x) = -A{x) . ^ ' ' 

Let us notice in passing that the "good" quark fields of definite helicity have 
a single nonzero complex component, which suggests that one can entirely 
get rid of the spinor indices. A general quasipartonic operator is built from 
the above parton fields 

0(^i, . . . ) = U^izi) n A{z,)l[i^l{zk)lli^l{zi) . (2.85) 

i j k I 

Its covariant analogue can be found expressing the (anti)holomorphic gauge 
potentials in terms of the field strength tensor, 

A = -l=a+i(F+^ - iF+^) = -^d^\F+y - iF+y) , 

A = ■j=d:^:\F+'' + = -^d^\F+y + iF+y) . (2.86) 

Contrary to "good" operators, the evolution of operators containing 
"bad" components of field operators does not have, in general, a simple 
pairwise pattern of the renormalization group mixing even at leading order 
in the coupling constant. A few exceptions do exist, however, for example 
in the case of the two-quark sector with Dirac structures 

i;{zin)il,j^^,^'[j')^{z2n). (2.87) 

After separation of twist-two and three components from these operators, 
one finds that, in the largc-Nc limit, both of them evolve autonomously with 
different pairwise Hamiltonians. The twist-two part is well known. The 
eigenvalues of the twist-three Hamiltonian coincide with the lowest anoma- 
lous dimensions for quark-antiquark-gluon operators introduced in Sect. 2.6 
( see [57-60]). 

The situation with more generic non-quasipartonic operators is an open 
issue. One of the major complications is that, in distinction with quasipar- 
tonic operators, the construction of the irreducible basis of non-quasipartonic 
operators is a nontrivial task. Up to now there exists no constructive pro- 
cedure for the classification of non-quasipartonic operators of arbitrary high 
twist. The representation theory used in the analysis of usual quasipartonic 
operators is not very handy since it does not allow one to get rid of the 
redundancy of the operators in the basis. The redundancy in the basis, in 
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turn, leads to complications in finding the eigensystem for the mixing ma- 
trix. Currently it can be solved only on a case-by-case basis. Finding a 
systematic procedure is an open question for further research. 

A particular case of certain non-quasipartonic local multi-gluon higher- 
dimension operators of the type 

L 

n^M.-.(0) (2.88) 
j=i 

in pure gluodynamics have been addressed recently and found to possess in- 
tegrable structures corresponding to the compact spin-one Heisenberg mag- 
net [61].^ Without loss of generality we can rephrase their analysis entirely 
in Euclidean space. The strength tensor can be decomposed into irreducible 
components 

F>^u = vtuF++vtuF- (2.89) 

with the help of 't Hooft symbols, 0(4) ~ SU{2) SU{2). The selfdual and 
anti-selfdual components transform as (1,0) and (0,1), respectively. The 
part of the RG Hamiltonian responsible for eigenvalues of autonomous com- 
ponents does not change the number of fields in the local gluonic operator 
By matching the coefficient of different irreducible components, ex- 



tracted by means of projectors P|^^ -^-^ for spin-j and parity P, to the avail- 
able one-loop calculations for gluonic operators up to dimension eight, the 
pairwise Hamiltonian was found to be [61] 

IHI12 = 7(F(2,o) +P(o,2))+F(i,o)+IP{o,i)-ll (Pj.o) +^P(o.o))+3IP{"i,i) • (^-^^^ 

The projection on the selfdual operators, i.e., built from products of F^, 
reduces the above Hamiltonian to 



If2 



7P(2,o)+F{i,o)-lllP{o,o), (2.91) 



where the projection operators extract maximal-spin, antisymmetric and 
trace components and have the following obvious representation 



\2,o)FM = \ {F^F^ + F^F^ - IS^^F^F^) , 

^(1,0)^ 



s pApB _ 1 jrAB pC j^C 



^ Other operators exist but do not affect the anomalous dimensions matrix as was demonstrated 
in [62] for L = 3, 4 by explicit calculation. 
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They can be easily related to the permutation PF^F^ = F^F^, trace 
KF^F^ = 5^^F^F^ and identity IF^F^ = F^Ff operators, 

Then, the pairwise Hamiltonian can be brought to the form 

ffi^ = 4Ji2 + 3Pi2 - 6Ki2 = 7 + 3si • S2 (1 - si • S2) , (2.92) 

where in the second equality we have used the representation in terms of 
spin-one SU{2) generators. This is a Hamiltonian of an exactly solvable spin- 
one Heiscnberg magnet which can be diagonalizcd by means of the Bcthe 
Ansatz [3,4,63]. It is interesting that the ground state of the corresponding 
spin-one chain is antiferromagnetic contrary to the ferromagnetic situation 
for the spin chains governing the renormalization of the scalar operators in 
J\f = 4 SYM theory. 

2.8. Conformal symmetry breaking 

Nothing is known at present about whether integrable structures survive in 
QCD beyond leading logarithms, in which case conformal symmetry is bro- 
ken by the trace anomaly of the energy-momentum tensor. The constraints 
imposed by conformal symmetry, and, in particular, by conformal Ward 
identities, appear to be nontrivial and can be used to advance perturbative 
calculations to higher orders. The idea has been to establish the following 
schematic structure of the perturbative series for a generic quantity Q 

Q = Q'^'^ + ^Maq , where AQ = power series in as . (2.93) 

Here Q'^"'^ is the result in the formal conformal limit, obtained by setting 
the /3-function to zero by hand. It has the full symmetry of a conformally- 
invariant theory. The extra term AQ can be perturbatively evaluated as 
a power series in ag and vanishes to leading order (LO). Note that the 
evaluation of the leading 0(a<j) contribution to AQ requires little effort, since 
it is sufficient to calculate the Nf proportional terms in /3o = 1 1 A^c/3 — 2A^j /3 
via quark bubble insertions. 

The possibility of the separation of the conformally-symmetric and /3- 
proportional contributions is by no means trivial. The reason is that al- 
though the dilatation anomaly to 0{as) accuracy receives contributions only 
from UV-divergent parts of relevant Feynman diagrams and is diagonal for 
conformal operators, the special conformal anomaly (to the same accuracy) 
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is affected by the finite parts of the same (one-loop) diagrams and is not di- 
agonal in a general, MS-like, renormalization scheme [64]. The nondiagonal 
entries in the special conformal anomaly can, however, be removed and the 
special conformal covariance (neglecting the /3-function) restored by perform- 
ing a finite renormalization [64]. This transformation defines the conformal 
subtraction scheme in which the separation in 1)2.93(1 becomes meaningful. 
Whether the remnants of integrability survive in this procedure, remains to 
be seen. 

3. Dilatation operators in supersymmetric gauge theories 

During the last year, a lot of new and interesting results have been obtained 
on the integrability properties of the = 4 SYM theory. These studies are 
important and provide hints on the long standing quest to find the corre- 
spondence between (supersymmetric) Yang-Mills theories at strong coupling 
and noncritical strings propagating on the curved background. Our aim in 
this section is to establish a connection between these new developments and 
the integrability of RG equations in QCD discussed in the previous section. 
In particular, the key question is to distinguish between phenomena that 
are already present in pure Yang-Mills theory (and only get enhanced by 
the supersymmetry), and ones that present genuinely new features of SYM 
gauge theories. 

In order to trace this connection, we introduce a superfield formalism 
in which all symmetries of the theory become manifest and calculations 
can be performed in a unified manner for different numbers of supercharges 
M = 0,1,2,4. The maximally-supersymmetric M = 4 SYM theory is a 
finite, four-dimensional conformal field theory [45,46,65,66], while the M = 
theory corresponds to pure gluodynamics. Unfortunately, no covariant 
superfield formulation for = 4 exists, which makes this case somewhat 
special. Because of this, we employ a noncovariant, light-cone superspace 
formalism due to Mandelstam [45] and Brink et al. [46] . This technique was 
used, e.g., to prove the finiteness of the = 4 SYM theory but is not very 
popular nowadays. Hence we begin this section by reviewing this beautiful 
framework. 

We will find it useful to introduce multi-particle single trace operators 
built from light-cone superfields at "super-light-cone" separations in the su- 
perspace Zk = where Zk are the usual "even" light-cone coordinates 
and 6^ are a set of Grassmanian "odd" coordinates. The one-loop dilatation 
operator for such operators can be found in terms of simple field displace- 
ments on the super-light-cone, similarly to the QCD case. Expanding such 
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operators in "even" and "odd" coordinates, we will be able to reproduce the 
integrability properties of QCD operators with multiple derivatives discussed 
in the previous section [10], on one side, and composite scalar operators, con- 
sidered in [25] on the other side, which are related to 5L(2,M) and 5*0(6) 
spin chains, respectively. Hence, these two results, that were found indepen- 
dently, appear to be just two sides of the same coin. We will then comment 
on the extension of these results to other operators. 

3.1. SYM theories in light- cone super space 

All super- Yang-Mills theories in four dimensions can be obtained through 
dimensional reduction from a higher-dimensional N = 1 gauge theory [67] 

S = j d'^x {- iF^^Ff ^ + i {d^5ab + gr'Af) V'"} , (3.1) 

describing the interaction of the D-dimensional gauge field A^'^ with the 
2[^/2]. component Majorana~Weyl fermion ^, both belonging to the adjoint 
representation of the gauge SU{Nc) group. In particular, the M = 2 SYM 
theory can be obtained from the D = 6 action, 

Cm=2 = tr |- iF^.F'^^ + 2iAiAaf P^A^ + 2 (P^0) (P^0) (3.2) 

- g^[<p,^f - V25eABA"^[0,Af] +^/2ge^^A^A[<A,Ai]| , 

(with A,B = 1,2), and the AA = 4 SYM theory follows from from D = W 
action, 

JOm=a = tr |- iF^.F'^'^ + i (P^0^^) {V^^ab) + Ig'i^'^'', 'P^^'Mab, ^cd] 
+ 2i~X^ACjfV>'X^ - v^gX'^^i^AB, Af ] + V2gX^A[<t>^'' , Ai]|, (3.3) 

(with A,B = 1, 4). Obviously, the first term on the r.h.s. of Eq. 1)3.2(1 gives 
the Lagrangian of a pure gauge, M = 0, theory, while the sum of the first 
two terms defines (for A = 1) the M = 1 SYM theory in the Weyl form. 
Thus, augmentation of supersymmetry increases the number of fundamental 
fields in the theory, which has profound consequences. The M = 0,1 SYM 
theories possess a non-vanishing beta-function which receives contributions 
at all orders of the perturbative expansion in the coupling. In = 2 theory 
the beta-function gets a contribution at one-loop only, while for AA = 4 
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the beta function vanishes to all orders. Hence the = 4 theory remains 
superconformally-invariant at the quantum level [45,46,65,66]. 

In order to construct a SYM theory on the light-cone, one starts with the 
component form of the action and fixes the light-cone gauge A^{x) = 0. In 
this gauge, one sacrifices four-dimensional Lorentz covariance in favor of a 
unified description of AA— extended SYM theories in superspace. The main 
steps in designing this formalism [45,46,68] consist of (i) making use of the 
decomposition of the elementary fields into propagating and non-propagating 
modes in light-cone time x+ ^, in the same way as indicated above in the 
context of QCD, (ii) integrating non-dynamical fields out in the functional 
integral, (iii) assembling the propagating fields into superfields, and finally 
(iv) casting the light-cone component action into the superspace form. 

The first two steps can be done for a generic D-dimensional M = 1 
Lagrangian without specifying the value of D. One finds that in the 
light-cone gauge A^{x) = 0, the action ()3.1|) involves two non-propagating 
fields [45,46,68]: one bosonic and one fermionic 

X^=r-- '^^^'r+^T>f^^;l , (3.4) 

which can be set to zero. Here X>^j = —dabd^ + ofacbA^ is the transverse 
covariant derivative. Solving 5"^ = = 0, one obtains the expressions for 
A°i and in terms of the propagating fields and substitutes them back 
into the SYM action. As the next step, one performs the reduction of the 
resulting action to four dimensions, giving rise to various SYM theories on 
the light-cone. Their explicit form for A/" > 1 is very lengthy and is not of 
interest per se, but rather as an intermediate result on the way to achieving 
the superspace form. 

We go directly to step (iii) and decompose all propagating, "physical" 
fields into definite helicity components. In the case of AA = 4 SYM,^ 
they include (anti)holomorphic, helicity-itl fields, A = {Ai + iA'i)l\pl and 
A = {Ai — iA2)/V2, built from two-dimensional transverse components of 
the gauge field, A±{x) = {Ai, A2), complex scalar fields of helicity and 
helicity ±1/2 components of Majorana-Weyl fermions, and A^. An im- 
portant property of the light-cone formalism, which makes it advantageous 
over the covariant one, is that the latter fields have only one non-vanishing 



''I.e., those that have or do not have the kinetic term with hght-cone time derivative 

' As we will see in a moment, the SYM light-cone theories with less supersymmetry can be obtained 

from the Af = 4 theory through the truncation procedure. 
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component = v^A^ and X\ = i \^Xa- As a consequence, one can de- 
scribe helicity-ibl/2 fermions by Grassmann- valued complex fields without 
any Lorentz index. Introducing the fermionic coordinates 9"^ possessing the 
helicity h9'^ = —^9^ and their conjugates 9a with Ma = ^^A, we can 
assemble the above fields into a single, complex, = 4 superfield [46] 

9^, 9a) = e^AS^ ^+ + 9^d^'~XAix) + 1^9^9''4>ab{x) 

-^eABCD^^^^^^A^(x)- ^eABCDe^e^9''9^d+A{x)y (3.5) 

which satisfies the chirality condition 

D'^'^ix, 9^, 9a) = {d-s^ - $(x, 9^, 9a) = . (3.6) 

It embraces all particle helicities, from —1 to 1 with half-integer step, and, 
therefore, (^{x,9^,9a) describes a CPT-self-conjugate supermultiplet. An- 
other advantage of the light-cone formalism is that the full supersymmetry 
algebra is reduced on the light-cone to a subalgebra which acts linearly on 
the "physical" fields and closes off-shell. In particular, in the AA = 4 SYM 
theory, the supersymmetric transformations on the light-cone read as 

M = C^Aa, (3.7) 
6<^^^ = ^ (^[^A^l - e^^^^icXn) , 

where [A,b]=ab-ba. 

Gauge theories on the light-cone with less or no supersymmetry can be 
deduced from the maximally supersymmetric AA = 4 theory by removing 
"unwanted" physical fields. In the superfield formulation this amounts to a 
truncation of the AA = 4 superfield, or equivalently, reduction of the number 
of fermionic directions in the superspace [69,70]. For instance, the M = 2 
chiral superfield is 

$(2) = $W|^3^^^^^ ^ e^^A^^^+i^d^^A{x)+9^d^^XA{x)+i9^9'^4){x)^ (3.8) 

with the identification = e^^(j) and A,B = 1,2. To get the M = I 
superfields one removes three odd coordinates 9^ = 9^ = 9^ = 0, and, finally 
for AA = all 0's in 1)3. 5() have to be set to zero. Notice that under this 
procedure the truncated M = 2, M = 1 and M = theories involve only half 
of the fields described by the AA— extended SYM theory and the other half of 
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the needed particle content arises from the complex conjugated superfields 
$ = For = 4 the two superfields ^ and ^> are related to each other 
through the reality condition [45,46] 

9^, 6 a) = -d:[:^D^D^D^D^^x, 9^, 9 a) , (3.9) 

but they are independent for < 2. 

The final step consists of constructing the action for the SYM theory in 
terms of the light-cone superfields. One starts with the = 4 theory, which 
involves a single superfield (|3.5|) . The action consists of the kinetic term 
as well as the cubic and quartic couplings $<^$<I>. Since light-cone su- 
persymmetry is realized linearly, each of the above three terms is separately 
invariant under the supersymmetric transformations (|H.7|) . Thus one can 
individually construct superspace invariants for each number of superfields 
and match them with the light-cone action in the component form to fix the 
overall normalization. This yields [46] 

- y^r''T'^^ (^d-\<^^d+^'')d~\¥d+^^) + I , (3.10) 

with $ given by (|3.9() . The SYM theories with less supersymmetry, M < 2, 
can be derived from 1)3. 1U() through the truncation procedure. For AA = 0, 1, 2 
the resulting light-cone action has the universal form [69] 

Sm<2 = jd^xd^9d^9!^axmdl--^'^'' (3.11) 

j\r(j\r+i) - 
where aj\f = (—1) 2 , Dj^ = Di . . . Dj^ and analogously for D . Here, in 

contradistinction with Eq. (|3.10|) . the superfields $ and ^ are independent 

on each other. 

3.2. Superconformal symmetry on the light- cone 

Our main objects of interest will be multiparticle single-trace operators built 
from light-cone superfields 

0(Zi,...,Zl) =tr{$(Zi)$(Z2)---$(ZL)}, (3.12) 
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where ^{Z) = ^"'{Z)t"' is a matrix valued superfield in the fundamental 
representation of the gauge group and Z denotes its position in the super- 
space with four even coordinates, x^, and 2Af odd coordinates, 9^ and 6a 
with A = l,...,J\f. In addition, we choose all superfields to be located 
along the light-cone direction in the four-dimensional Minkowski space de- 
fined by the light-like vector (with = 0) defined as in Eq. (|2.1|) . so 
that n ■ A = Aj^ = 0, cf. H2.2|) . Similarly to the QCD case, the posi- 
tions of the superfields on the light-cone are parameterized by real numbers 
= zn^. One can further simplify the notations by taking into account 
that a ^—dependence of the chiral superfields can be absorbed into the re- 
definition of z— variables 

9^, 9a) = ^{{z + \9- 9)n^, 9^, O) . (3.13) 

Redefining z + ^9 ■ 9 ^ z, we shall specify the position of the superfield in 
the superspace as Z^ = {zk,9^) so that ^{Z^) = <I>(2fcn^, 0). 

The single-trace operators (|3.12|) represent a natural generalization of 
nonlocal light-ray operators in QCD, cf. Eq. ()2.6|) . To obtain the latter it is 
sufficient to expand 0(Zi, . . . , Zl) in powers of odd variables 9^'... . As 
in QCD, nonlocal operators ()3.12|) serve as generating functions for Wilson 
operators with the maximal Lorentz spin and minimal twist equal to the 
number of constituent fields L. A question arises as to whether this set 
of Wilson operators is complete. In = 4 SYM theory there is only one 
independent chiral superfield ^{Z) and, as a consequence, the operators 
H3.12|) generate all Wilson operators of twist— L built from L fundamental 
fields. For M < 2, the superfields ^{Z) and ^{Z) are independent of each 
other and, in addition to the operators in ()3.12() . one can introduce "mixed" 
operators built from both superfields. This means that in the M = 0,1 
and 2 SYM theories, the operators (|3.12|) only generate a certain subset 
of the existing Wilson operators. The reason why we will concentrate on 
this subset is that it inherits integrability properties found in QCD in the 
sector of maximal helicity operators. Namely, as we shall argue below in this 
section, the one-loop dilatation operator acting on the space of single-trace 
operators (|3.12jl can be mapped in the multicolor limit into a Hamiltonian 
of a completely integrable Heisenberg spin magnet with the symmetry group 
being a collinear SL(2\J\f) subgroup of the full superconformal group [24]. 

To explain this correspondence, we would like to recall the symmetries of 
the SYM theories. At the classical level, a four-dimensional supersymmetric 
Yang-Mills theory with M supercharges is invariant under the SU{2,2\J\f) 
superconformal transformations [67]. The corresponding algebra contains 
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both bosonic (even) and fermionic (odd) generators. The even generators are 
the fifteen generators of the 50(2,4) conformal algebra famihar from QCD 
plus the chiral R charge and the SU{J\f) internal "flavor" symmetry gener- 
ators Ta^ (for Af > 2 only). The odd generators are the generators of su- 
persymmetric and superconformal transformations, QaA, and S;^, S^, 
respectively. Since the superfields entering 1)3.12(1 "live" on the light-cone, 
we expect, by analogy with the QCD case, that the full superconformal 
group will be reduced to its collinear subgroup. We know from Sect. 2.1 
that in Yang-Mills theories with no super symmetry, the collinear subgroup 
is generated by three bosonic operators L^'*^ which form the SL{2, M) alge- 
bra. Supersymmetry enlarges the 5L(2,M) subgroup by adding the follow- 
ing generators: the chiral charge R, SU{J\f) charges T^^, transverse boosts 
Mi2 and, finally, the "good" components of the fermionic charges Qa, Q^, 
and Sa- Being combined together, these operators form the SL{2\J\f) 
algebra [24]. 

To one-loop accuracy, the generators of the SL{2\J\f) subgroup can be 
realized as linear differential operators, [G, = G{dz,dg)^, acting on the 
light-cone coordinates of the chiral superfield ^(Z) [24,71]. The correspond- 
ing expressions for the generators of the SL{2; M) collinear subgroup gener- 
alize similar expressions in the QCD case, Eq. (|2.18|) . 

L- = -d, , L+ = 2jz + z^d, + z{9-de) , L^" = j + zd, + 1 (0 • 9^) , 

(3.14) 

and the twist operator is E = t. The four odd generators can be represented 
by the following differential operators 

W^'- = e^d,, W^'^ = 9^[2j + zd, + {9 ■ de)] , VX=deA, VX = zd^A . 

(3.15) 

Finally, the S\J (M) generators and linear combination of the abelian 
i?— charge and transverse boosts M12 are realized as 

TB^ = e^deB-^6^{e.de) , B = -j - {l - ^) {d ■ dg) . (3.16) 

In these expressions, dz = d/dz and 9 ■ dg = 9'^d/d9^. The twist t and 
conformal spin j are the parameters defining the representation. For the 
chiral superfield ^{Z) entering the Lagrangian of AT-extended SYM theories, 
Eqs. (IXTU)) and (|57n)) . one has j = -\ and t = 1. 

The hght-cone operator 0(Zi, . . . , Zl), Eq. (|3.12|) . serves as the generat- 
ing function of an infinite tower of Wilson operators with an arbitrary num- 
ber of derivatives, when expanded in the bosonic separations, and of different 
field content, when expanded in fermionic coordinates. The most efficient 
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way to renormalize an infinite set of local Wilson operators is to study the 
renormalization properties of the light-ray operators 0{Zi, . . . , Z^) them- 
selves. Superconformal invariance imposes severe restrictions on the possible 
form of the dilatation operator in the AA— extended SYM theory acting on 
the space of light-cone operators H3.12() . For the operator 0{Zi, . . . , Z^), 
the generators of the SL{2\J\f) coUinear subgroup are given by the sum of 
single-particle generators, Eqs. (|H.14j) - acting on the coordinates in 

the superspace Z^ = {zk,0^) with k = 1, . . . , L. In particular, the total twist 
of the operator 0{Zi, . . . , Zl) is equal to the number of superfields L. It is 
well known from QCD [12] that to one-loop accuracy the RG evolution of 
local quasipartonic operators preserves the number of particles. In addition, 
the mixing of single trace operators with operators involving larger number 
of traces is suppressed by powers of Taken together, these two ob- 

servations imply that to one-loop accuracy the set of light-cone operators 
0{Zi, . . . , Zl) is closed under renormalization in the multicolor limit. The 
corresponding one-loop dilatation operator can be represented in the light- 
cone superfield formalism as a quantum mechanical Hamiltonian EI leading 
to the following RG equation 

( d d ^ a^N 

|m ^ + + L7(<7)| 0^(Zi, . . . , Zi) = [M • 0«] (Zi, . . . , Zz.). 

(3.17) 

The renormalized operator generates finite Green functions and is ex- 
pressed in terms of the bare operator ()3.12|1 with an infinite multiplicative 
counterterm. Here 7(5) is the anomalous dimension of the superfields in 
the light-like gauge Aj^{x) = 0. Thanks to the Ward identity it is related 
to the beta-function of the theory ^{g) = [3{g)/g. The Hamiltonian H in 
Eq. (|3.17j) is an integral operator acting on the coordinates of the superfields 
in the superspace. Superconformal invariance implies that the dilatation 
operator IH commutes with the generators of the SL{2\N) supergroup in 
Eqs. (nmjl - dTTHll . 

3.3. One-loop dilatation operator 

To one-loop order, the dilatation operator H has a pairwise structure. In 
addition, only the nearest-neighbor interactions survive in the multicolor 
limit, leading to 

EI = IHi2 + IHI23 + ---+Hli. (3.18) 

The two-particle interaction kernel ^k,k+i receives contributions from the 
Feynman diagrams shown in Fig. 0J They involve cubic and quartic inter- 
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Figure 4. Feynman diagrams contributing to tlie one-loop dilatation operator for the two chiral 
superfields on the light-cone. 



action vertices which are uniquely defined by the hght-cone actions (|3.1U|) 
and ()3.12|) . As a consequence of superconformal symmetry the kernel M^^k+i 
has to commute with the generators of the SL{2\J\f) superalgebra acting on 
the superfields with labels k and A; + 1 and, therefore, it is a function of the 
two-particle quadratic Casimir operator 



^k,k+l 



(3.19) 



where a general expression for the quadratic Casimir in the SL{2\N) algebra 
is [72] 



+ {N- + 



(3.20) 



M-2 



1 rpB rpA 



For = 0, this expression coincides with Eq. 1)2. 21() . In that case, the SL{2) 
invariance allows one to obtain the general form of the two-particle kernels 
in QCD, Eq. ^(T^ . 

We remind the reader that in the light-cone superspace formalism, the 
Yang-Mills theory without supersymmetry can be obtained from the one 
with M supercharges by putting the odd coordinates to zero, 0^ = 0. For 
this reason the integral operator 1)2. 3U() can be considered as the special limit 
of a general two-particle kernel M.k^k+i for 9^ = O^j^^ = 0. Making use of 
the odd generators in the SL{2\M) algebra one can "lift" the operator (|2.3U|) 
into the superspace and restore the dependence on the odd coordinates. In 
this way, one finds a general form of the two-particle SL{2\N) invariant 
evolution kernel 



[Mi2-0](Zi,Z2,...,Zl) (3.21) 

da Hapf^'^ 0{aZi + aZa, +PZ2,..., Zl) , 

10 Jo 



where j = —1/2 is the conformal spin of the chiral superfield ^{Z) entering 
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into the light-cone operator ()3.12|) . Here the kernel M12 acts on the coordi- 
nates of the superfields ^{Zi) and ^{Z2) and displaces them simultaneously 
along bosonic and fermionic directions in the superspace Z = {z,9^). 

The evolution kernel ()3.21|) involves an arbitrary function if of the har- 
monic ratios. Its form is not fixed by superconformal invariance and is 
determined by the dynamics of the underlying SYM theory. The one-loop 
calculation of super-Feynman graphs, displayed in Fig. |1J gives (/?(^) ~ 
This suggests that the two-particle kernel looks like [24] 



This relation admits an obvious generalization to Yk,k+i- However, 
Eq. (|3.22|) cannot be identified as the two-particle dilatation operator since 
the integral is divergent for a ^ 0. The divergence arises due to a negative 
value of conformal spin j = —1/2 of the light-cone superfield - see Eqs. (|3.5() 
and ()3.8() - which can be traced back to the fact that the lowest component 
of the superfield is a nonlocal operator d^^A{z). This means that expanding 
the nonlocal operator 0(Zi, Z^) around the origin Zi = ... = Zl = 0, in 
terms of local operators one obtains both Wilson operators and "spurious" 
operators involving nonlocal fields d^^A{0), ^+(0) and d^^X{0). To elimi- 
nate the latter it is sufficient to subtract the two first terms in the expansion 
of the superfield ^{Z) around Z = 0: 



where the superscript indicates that the superfield $^(Z) does not generate 
nonlocal fields. Let us define the operator H that projects the nonlocal 
operator ()3.12|) onto the subspace of Wilson operators 



0^(Zi,...,Zl) =nO(Zi,...,Zi) =tr{cI>*(Zi)---cDW(ZL)}, (3.24) 



so that n2 = H. Since Wilson operators mix among themselves and cannot 
mix with spurious operators (but the opposite is possible!), one expects 
that nonlocal operators 0^(Zi, Zl) evolve autonomously under the RG 
transformations. Indeed, decomposing the operators as 



0(Zi, Zl) = 0^(Zi, Zl) + (1 - n)0=''"^(Zi, Zl) (3.25) 



and substituting this relation into (|3.22l) . one finds that the operators 
0^{Zi, Zl) obey the same relation (|3.22|) but with the two-particle kernel 




0(aZi + aZ2, Z2, Zl) - 0(Zi, aZa + aZi, Zl)} . 



<^{z) <^^{z) = <^{z) - $(0) - z ■ az^>(o) 



(3.23) 
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replaced by [24] 

Vi2 -^mY2 = nvi2n . (3.26) 

It is straightforward to verify that the integral operator is well defined, 
that is, the above-mentioned divergences for a ^ do not appear in the 
expression for ]HI]|^0(Zi, Zl). Thus, as far as the renormalization of Wil- 
son operators is concerned, the one-loop dilatation operator is given in the 
multicolor limit by Eq. ()3.18|) with the two-particle kernel IH[^^_(_j^ defined in 

The expression obtained for the dilatation operator does not exhibit an 
explicit dependence on M. This dependence enters entirely through the 
number of odd directions in the superspace, Z = (z,6'^). This means that 
for nonlocal light-cone operators ()3.12|) the dilatation operator has the same, 
universal form inAA = 0,AA = 1, M = 2 and M = 4 SYM theories on the 
light-cone. The important difference between these theories is that for = 4 
the operators (|3.12() generate all Wilson operators of the twist L, whereas for 
M < 2 they capture only a subsector of Wilson operators. To describe the 
remaining operators, one has to consider single-trace operators built from 
both chiral and anti-chiral superfields, $ and respectively. Remarkably 
enough, the dilatation operator for such mixed nonlocal operators can be 
deduced from the dilatation operator (|3.22|) in maximally supersymmetric 

= 4 theory via the truncation procedure described above. 

3.4. Dilatation operator as a SL{2\Af) Heisenberg spin chain 

Let us consider the properties of the one-loop dilatation operator 1)3. 22(1 
in more detail and demonstrate how one can use the obtained expres- 
sions to evaluate the anomalous dimensions of various Wilson operators in 
AA— extended SYM theories. 

To begin with, consider 1)3.22(1 for AA = 0. In this case, the nonlocal 
operator (|3.12() is reduced to the product of holomorphic components of 
gauge fields d^^A{z). To avoid spurious operators it is convenient to go over 
to the nonlocal light-cone operators 

Og{z,,. ■■,zl) = {d,, . . . d,^fO{zi, ...,zl) =ti {d+A{zi) ■ ■ ■ d+A{zL)] . 

(3.27) 

Remember that we have adopted the light-like axial gauge Aj^{x) = 0. 
In the covariant formulation the field dj^A{z) is replaced by n^{F^±(z) — 
iFf^±{z)). Recall that the field A{z) carries the helicity —1 and, therefore, 
Og{zi, . . . , zl) is the generating function for pure gluonic Wilson operators 
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of maximal helicity —L. 

Substituting Eq. (1^771) in one gets 

[ml^-03]izi,Z2,...,ZL)= —\20'^{azi + az2,Z2,...,ZL) (3.28) 

Jo a 

-a'^0^{azi + QZ2, Z2,. ■ ■ , zl) - oP'O^ {zi,az2 + azi, . . . , zl)| . 

This result should be compared with the similar expression for the dilatation 
operator for quarks with the aligned helicities, Eq. (|2.26|) . The only difference 
is in the power of a in front of the second and third terms. In general, this 
power is determined by the conformal spin of the field, ~ q^-'"^. For the 
quark field, one has jq = 1, whereas for the gauge field strength, jg = 3/2, 
and for scalar fields, js = 1/2. 

To diagonalize the dilatation operator, one has to repeat the same steps 
as described in Sect. 2.2 in the case of aligned helicity quark operators. 
Namely, one substitutes 

(Q3{zi,Z2,...,Zl)^{zi-Z2T = Z^2 

into (|3.28|1 and obtains 

Mfa 2^2 = 2[V'(n + 3) - ^(1)1^^2 • (3-29) 

This relation defines the anomalous dimension of the conformal operator 
built from two gauge fields of the same helicity 

OgiQ) = d+A(Q){id+Td+A(Q) + . . . . (3.30) 

Here ellipses stand for terms involving total derivatives. The conformal 
iS'L(2;M) invariance fixes them uniquely and allows one to write down a 
closed expression in terms of Gegenbauer polynomials [15]. The conformal 
spin of the operator ()3.3U|) equals J12 = n + 2jg = n + 3 and, therefore, the 
two-particle dilatation operator 1)3. 28(1 can be written in the SL{2) invariant 
form 

IHf2 = 2[V'(Ji2)-^(l)] (3.31) 

with J12 being the two-particle conformal spin. In contrast with the quark 
operator, the corresponding S'L(2;M) representation has jg = 3/2. 

Combining together ()3.31|) and 1)3. 18() . we conclude that the dilatation 
operator in the sector of maximal helicity gluon operators can be identi- 
fied as the Hamiltonian of the S'L(2;M) Heisenberg magnet. The length of 
the spin chain is equal to L and the value of the spin in each site is given 
by the conformal spin of the gauge strength tensor jg = 3/2. Notice that 
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the maximal-helicity gluonic operators evolve autonomously and the corre- 
sponding one-loop mixing matrix is sensitive to the gluonic sector of the 
gauge theory only. This implies that the one-loop dilatation operator in the 
sector of maximal-helicity gluonic operators is the same in all AA— extended 
SYM theories. 

The Wilson operators built from six scalars (l)^^{x) with no derivatives 
provide another limiting case which reflects the specifics of the = 4 SYM 
theory. It proves convenient to introduce their linear combinations [67] 

4>^(x) = ^^f^AB{x) (3.32) 

where j = 1, ... ,6 and "^f^ are the chiral blocks of Dirac matrices in six- 
dimensional Euclidean space. Then, the composite scalar operators look like 



Onn...j, = tr{(/>,,(0)</.,,(0) . . . cPj^m . 



(3.33) 



To one-loop order, these operators evolve autonomously under the RG trans- 
formations and the corresponding dilatation operator has the form 1)3. 18() 
with the two-particle kernel given by a finite-dimensional matrix 



(3.34) 



The scalar field can be obtained from the N" = 4 superfield p.5|) by means 
of a simple projector 



9^=0 



(3.35) 



Using this projection, the mixing matrix Vj_^j^ can be extracted from the 
general expression for the = 4 one-dilatation operator given in Eqs. 1)3. 26() 
and 



31J2 



{Zl,Z2 



Zi=Z2=0 



The calculation of fermionic derivatives then yields 



iii2 



32 



(3.36) 



(3.37) 



This matrix acts on the isotopic indices of two neighboring scalar fields in 
H3.33() and has the size 6^ x 6^. Hence, the one- loop dilatation operator for the 
scalar composite operators 1)3. 33() is given by the 6^ x 6^— matrix. As in the 
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previous case it has a hidden symmetry - this matrix can be mapped into 
the Hamiltonian of the completely integrable Heisenberg S0{6) ~ S'C/(4) 
spin chain of length L [25] . 

We see that in the Af = 4 SYM theory, the integrability phenomenon 
occurs in two different sectors of Wilson operators. In both cases, one en- 
counters the Heisenberg spin chain of length L but with different symmetry 
groups. For the maximal helicity Wilson operators with arbitrary number of 
derivatives, the SL{2;M) symmetry group arises and it is related to confor- 
mal symmetry on the light-cone. For scalar operators with no derivatives, 
the arising SU{4:) symmetry group is related to the i?— symmetry; that is, 
rotations in the isotopic space. From the point of view of the superspace, 
Z = {z,6^), the two groups act separately on the bosonic and fermionic 
coordinates of the superfields and as such they are subgroups of a bigger, 
SL{2\M) group. Since the SL{2; M) and SO{Q) spin chains emerge as special 
limits of a general expression for the one-loop dilatation operator Eq. 1)3.22(1 . 
a question arises about integrability properties of the latter operator as a 
whole. Quite naturally, the one-loop dilatation operator can be mapped into 
a Hamiltonian of the Heisenberg spin magnet with the SL{2\J\f) symmetry 
group. The latter is defined in ()3.18|) with the two-particle kernels given by 
the familiar expression [24, 71] 

11^2 = 2 [^(Ji2) - V(l)] , Jl2{Jl2 - 1) = ^?2 (3.38) 

in which J12 is the two-particle SL{2\J\f) spin and the quadratic Casimir was 
introduced in ((3.2U|) . 

For = 4 this (super) spin chain describes renormalization of arbitrary 
Wilson operators in the underlying SYM theory built from "good" compo- 
nents of the fundamental fields. For AA < 2, it covers only a special subsector 
of Wilson operators built from fundamental fields entering into the expres- 
sion for the chiral superfield, Eq. (|3.8|) . In particular, at M = this sector 
corresponds to maximal helicity gluonic operators. For the remaining Wilson 
operators generated by mixed products of chiral and antichiral superfields, 
the dilatation operator does not exhibit integrability properties. As was il- 
lustrated in Sect. 2.5 in the case of non-maximal helicity quark operators, 
the additional terms in the dilatation operator responsible for breaking the 
integrability generate a mass gap in the spectrum of the anomalous dimen- 
sions. 

The findings of Ref. [23] suggest that in the AA = 4 theory, not only 
quasi-partonic operators, addressed above, but all operators obey integrable 
renormalization group equations. Thus the spin chain discussed above gets 
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enhanced to the SU{2,2\4) superchain. As we have emphasized in Sect. 12.71 
this is a very nontrivial statement in the QCD context. It suggests that 
operators of different twists are intertwined in a very subtle manner. The 
superconformal symmetry does indeed relate operators of different twists. 
However, this is not enough to conclude that the anomalous dimensions of 
non-quasipartonic operators are expressed via the quasipartonic operators, 
since this is not the case in the most straightforward supersymmmetric gen- 
eralization of QCD, the Af = 1 SYM. The presence of the elementary scalar 
fields in the = 4 model does matter, so that multiple full^ supersymmetric 
transformations of the operator (p"^ (p"^ , having the integrable form of the 
anomalous dimension (|3.29j) will yield combinations of operators with differ- 
ent twists but having the same anomalous dimension. This is not the case in 
QCD where anomalous dimensions of different twist operators are believed 
to be dynamically independent. This issue deserves thorough further study. 

There have been several recent related developments that we cannot ad- 
dress here in detail. In particular, the anomalous dimension matrices in 
different autonomous sectors in = 4 SYM were studied in Refs. [73-76], 
mostly in the different subsectors of scalar operators. Another very impor- 
tant question is whether the integrability of the dilatation operator persists 
at higher orders in the loop expansion. Several positive demonstrations that 
this does happen have been given in certain scalar subsectors [77-79] . To ex- 
tend these findings to the full dilatation operator, the analysis of the SL{2) 
subsector has to be performed. The two-particle part of the corresponding 
dilatation operator was calculated in Refs. [80] up to three-loop order. 

Ultimately, one has to understand what symmetry of the quantum field 
theory is responsible for the appearance of integrability [81, 82] and what 
part of it survives after decreasing the amount of supersymmetry [83-85], as 
we explained for one particular case above. 

4. High Energy Scattering 

In the previous section, we have considered the perturbative QCD dynam- 
ics of partons constrained to small transverse separations, which can be 
achieved by considering "hard" scattering processes in generalized Bjorken 
kinematics. When the transverse distance becomes smaller than any other 
scale in the process, the corresponding dependence can be treated by 
renormalization-group methods and disappears from the dynamics. The 
problem becomes essentially one dimensional: the relevant degrees of free- 



' I.e., not restricted to the light-cone as in Eq. 13.71 . 
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dom correspond to partons "living" on the light-ray and the interaction 
between them amounts to displacement along the light-ray, with a certain 
weight-function. We were able to analyze this interaction using a Hamil- 
tonian formulation and found its hidden integrability property. The scale 
dependence in this case corresponds to "time" evolution in the quantum- 
mechanical picture. 

The present section is devoted to perturbative QCD dynamics in the 
high-energy limit. In this case, the invariant energy s = {pA + Pb)^ of 
colliding hadrons A and B is considered to be the largest scale and, in 
the limit s — > oo, the energy dependence corresponds to a renormalization 
group flow for a dynamical system "living" in the two transverse dimensions 
orthogonal to the scattering plane. The relevant degrees of freedom corre- 
spond to reggeized gluons (see below) with assigned transverse coordinates, 
and the interaction between them leads to a displacement in the transverse 
plane. Compared to the "hard" scattering situation, the longitudinal and 
transverse degrees of freedom change places, and the dynamics becomes two 
dimensional in the transverse plane instead of being one dimensional in a 
given light-ray direction as in the former case. In what follows, we will sketch 
the Hamiltonian approach to high-energy scattering and establish its hidden 
symmetry properties. The main result of this analysis is that the system of 
interacting reggeized gluons can be mapped into a quantum XXX Heisen- 
berg magnet with the SL(2,C) symmetry group [7-9]. This means thatthe 
high-energy asymptotics of scattering amplitudes in QCD in a certain limit 
can be calculated using the Quantum Inverse Scattering Method [32] . 

4.1. The BFKL Hamiltonian 

The high-energy asymptotics of the scattering amplitude A{s, t) is studied 
most conveniently through the Mellin representation which corresponds to 
an expansion in partial waves with complex angular momentum 

A{s,t) = is —s'^A{u,t). (4.1) 

Here, the integration contour goes to the right of all singularities of v4(a;, t) on 
the complex w— plane. The high-energy asymptotics of A{s, t) is determined 
by the singularities of the partial wave amplitudes: If ^(w, t) ~ 1/{uj — ujq (t) ) 
then A{s,t) Poles in the a;-planc arc called rcggeons, and the 

position of the pole as a function of the transfered momentum, LOo{t), is 
called the reggeon trajectory. Prom the point of view of QCD perturbation 
theory, finding the high-energy asymptotics of the amplitude corresponds 
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to the resummation of logarithms of the energy ~ (aglns)™. By virtue 
of the Mehin transformation (|4.H) . the expansion over (oslns)™" is traded 
for the expansion over (a^/a;)™' with coefficients that depend on t, and the 
power-hke energy-dependence arises because the resummed expression for 
A{as/uJ,t) develops a nontrivial singularity for positive uj = iOo{t) in the 
sum to all orders in as- 

The partial wave amplitude A{uj, t) can be written in the so-called impact- 
parameter representation 

Aiuj,t) = j ShoS'^^''^ j (fhd%' ^A{h-hQ)T^{b,b')^B0) 

= I dXe'^'^'^H^bomm), (4.2) 

where the impact factors ^A(b) and stand for the parton (gluon) dis- 

tributions as functions of transverse coordinates b = {bi,b2, ■ ■ ■ ,bn} and 
b' = {b[,b2, ■ ■ -b'^} in the colliding hadrons A and B, respectively, and 
Tuj{b, b') is the scattering (partial wave) amplitude for a given parton config- 
uration. The existence of such a representation follows from the space-time 
picture of the high-energy collisions: The two Lorentz-contr acted hadron 
discs approach each other at the speed of light and the interaction between 
them takes a finite time in the center-of-mass frame. The preparation of the 
initial states is separated in time from the interaction between hadrons and, 
therefore, the transverse coordinates of parton constituents in hadrons can 
be considered as "frozen" during the interaction. Note that the factorized 
form in (|4.2I) implies that the structure of singularities in the complex uj- 
plane does not depend on the parton distribution in hadrons, but rather on 
general properties of the gluon interaction in the i-channel. 

A nontrivial analysis [6] shows that the interaction effects modify the 
propagators of the i— channel gluons in such a way that they develop their 
own Regge trajectory; such t— channel gluons "dressed" by virtual correc- 
tions are called reggeized gluons. The reggeized gluons with given transverse 
coordinates provide one with the relevant degrees of freedom for the study of 
high energy interactions and play a role similar to that of the fundamental 
fields with given light-ray coordinates for light-cone dominated processes. In 
the course of interaction with each other, the reggeized gluons change their 
two-dimensional transverse coordinates and the color charge. The partial 
waves T^(6, b') can be classified according to the number of reggeized gluons 
propagating in the t-channel; the minimum number of two gluons is required 
in order to get a colorless exchange. Moreover, it can be shown that retain- 
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ing only two t-channel reggeized gluons corresponds to the resummation of 
energy logarithms to the leading logarithmic accuracy. In this approxima- 
tion the amplitude Tt^(6i, 62; ^'ii ^2) depends on four transverse coordinates 
and satisfies a Bethe-Salpeter-like equation - the so-called BFKL equation 
- in the operator form 

LoT^ = ¥(,0) + ^ Mbfkl , (4.3) 

TT 

where T^*^^ corresponds to the free exchange of two gluons. The solution of 
this equation can formally be written as 

= (^^ - ^ Mbfkl) ' t£o) , (4.4) 

whence one concludes that the singularities of in the w-plane are deter- 
mined by the eigenvalues of the BFKL operator 

PBFKL • ^a] {blM = Ea ^aihM) , (4.5) 

with a enumerating the solutions. The high-energy behavior of the scatter- 
ing amplitude is governed by the right-most singularity of T^, which cor- 
responds to the maximal eigenvalue maXoE^Q. Equation (|4.5|) has the form 
of a Schrodinger equation for a system of two interacting particles on the 
two-dimensional plane. Such particles can be identified as reggeized gluons 
and the eigenstates ^q(6i, 62) have the meaning of the wave functions of the 
color-singlet compound states built from two reggeized gluons. 

The BFKL operator Hbfkl has a number of remarkable properties which 
allow one to solve the Schrodinger equation 1)4. 5(1 exactly [86,87]. First of 
all, Hbfkl splits into the sum of two operators acting on the holomorphic 
and the antiholomorphic coordinates 

bj = {xj,yj} =^ Zj = Xj + iyj , zj = xj - iyj , (4.6) 

(with j = 1,2) on the transverse plane: 

Hbfkl = ^2+^2 (4.7) 

with 

n2 = ln(zi2) d,, + ln(zi2) d,, + HdzM - 2^{1) , (4.8) 

where Z12 = zi — Z2 and 7^2 is given by a similar expression in the 2— sector. 

Another remarkable property of Hbfkl is that it is invariant under the 
conformal SL{2; C) transformations of the reggeon coordinates on the plane 
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^^^azk + h^ (ad-6c=l). (4.9) 

czk + a 

Consider the generators of the holomorphic SL{2, C) transformations 

Lk - = -dzk , Lkfl = Zkdzk , Lk,+ = zldz^ , (4.10) 

and the corresponding antihofomorphic generators L^^-, L^fl and Lk^j^ given 
by similar expressions with Zj. replaced by z^, with A; = 1,2 enumerating 
particles. By inspection one finds that Hbfkl commutes with all two-particle 
generators 

[Hbfkl, ^i,a + L2A = Pbfkl, U^a + U^a] = (4.11) 

with a = +, — , 0. This implies that Hbfkl only depends on the two-particle 
Casimir operators of the SL{2,C) group 

Ll2 = -{zi-Z2fd,,d,,, Ll2 = -{zi-Z2fd,,d,,, (4.12) 

and leads to 7^2 = "^2 (-^12) Siud TI2 = 7^2 (-^12)- & consequence, solu- 
tions of the Schrodinger equation 1)4. 5|) have to be eigenstates of the Casimir 
operators 

Ll^^n,u = hih-l)^n,u, Pu^n,u = h(h-l)^n,u- (4.13) 

Here a pair of complex conformal spins is introduced 

1 "t" ^ — 1 — Tl 

h = — h iv and h = — h iv (4-14) 

with a non-negative integer n and real u that specify the irreducible (princi- 
pal series) representation of the 5L(2,C) group to which 'ifn,u belongs. The 
solutions to Eqs. 1)4. 13(1 read as 

(l+ra)/2+ii^ / - \ (l-n)/2+ii^ 

^) , (4.15) 

2^10 ^;2o/ 



\Zl0Z20j V 



where Zjk = Zj — z^ and bo = {zo,zo) is the collective coordinate, refiecting 
the invariance of Hbfkl under translations. We recall that the solution 
1)4. 15() defines the wave function of the color-singlet compound state, which 
is built from two reggeized gluons with the coordinates bi = {zi,zi) and 
&2 = (-225^2)- The integer n fixes the two-dimensional Lorentz spin of the 
state, the real valued u gives the scaling dimension £ = 1 + 2iv, and the 
two-dimensional vector 60 sets up the center-of-mass coordinate of the state. 
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To find the eigenvalues in Eq. (|4.5I) . one substitutes the wave function 
(|4.15|) into the Schrodinger equation (|4.5|) and uses the exphcit form of the 
BFKL kernel leading to [6] 

En,u = 20(1) - f ^ + f ^ - iu) . (4.16) 



Its maximal value, max En^u = 4 In 2, corresponds to n = z/ = 0, or equiva- 
lently h = h = 1/2. It defines the position of the right-most singularity of 
the partial wave amplitude and is translated into the asymptotic behavior 
of the scattering amplitude in the leading logarithmic approximation, 

A{s,t)^zs'+'^^'-\ (4.17) 

known as the BFKL pomeron.'^ 

Using ()4.16|1 one can reconstruct the operator form of the BFKL kernel 
Hbfkl on the representation space of the principal series of the SL{2, C) 
group 

Hbfkl = ^ [H{Ji2) + H{Ji2)] , H{j) = 2tP{l)-i;{j)-i^{l-j) , (4.18) 

where, as before, the two-particle spins are defined as = «/i2(«^i2 — 1) 
and = Ji2{Ji2 — !)• Notice that we already encountered the similar 
Hamiltonian in Sect. 2 and found that it gives rise to complete integr ability 
for the three-particle evolution equations for the helicity— 3/2 baryon distri- 
bution amplitudes. It turns out that the BFKL kernel 1)4. 18|) has the same 
hidden integrability. In order to appreciate this, we have to consider states 
containing more than two particles. 



4.2. Multireggeon compound states 

In the leading logarithmic approximation, the scattering amplitude 1)4. 17(1 
grows as a power of the energy and violates unitarity constraints. To re- 
store unitarity, subleading corrections must be included. Going beyond the 
leading logarithmic approximation, one has to take into account contribu- 
tions to the scattering amplitude in Eq. (|4.2|) with more than two reggeized 
gluons propagating in the t-channel. One can argue [88, 89] that, in or- 
der to restore the unitarity of the scattering amplitude in the direct {s— 
and t—) channels, it suffices to retain contributions where the number of 

''Due to accumulation of the energy levels -Bn,i/ around n = v = 0, the corresponding 
Regge singularity is not a pole but a square-root cut so that 14.171 is modified by additional 
(oa Ins)"-*-/^— factor. 
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the t-channel gluons is conserved, i.e., neglect the three-reggeon vertices in 
the Regge theory language. With this assumption, the amplitude AN{s,t) 
corresponding to the exchange by reggeized gluons satisfies the Bartels- 
Kwiecinski-Praszalowicz (BKP) equation [88,90], which is the generalization 
of the Bethe-Salpeter equation (|4.3() to the A^— gluon scattering amplitude. 
Singularities of in the complex w-plane are determined by the eigenvalues 
of the BKP operator Hat. They can be thought of as energies of the A^-gluon 
compound states 

El7V^iv,M(&i, ■■■M) = Sjv,m^7V,m(&i, ■■■M), (4.19) 

with {q} being some set of quantum numbers parameterizing the solutions, 
and bn the two-dimensional transverse positions of n— th reggeized gluon. 

In the large Nc limit, the relevant Feynman diagrams have the topology 
of a cylinder and Hat reduces to the sum of terms corresponding to pairwise 
nearest-neighbor BFKL interactions: 

^n = IJ:mZI\ (4.20) 

k=l 

with periodic boundary conditions H^^^^j^ = H^^^^^. The holomorphic 
separability of Hbfkl implies the same property for 

MN = nN + nN (4.21) 

with the operators TCn and TYat acting on z— and z— coordinates, respec- 
tively. The conformal symmetry is also inherited by them, so that each of 
the two-particle Hamiltonians can be written in terms of the corresponding 
two-particle Casimir operators 

1 ^ 

= 2 ^ [^^(^fc,fc+i) + HiJk,k+i)\ , (4.22) 

k=l 

where H{j) is given in terms of -^-functions (see Eq. 1)4. 18() ^. 

As a consequence, the (2+l)-dimensional Schrodinger equation 1)4. 19(1 
can be replaced by the system of holomorphic and antiholomorphic (1+1)- 
dimensional Schrodinger equations for the Hamiltonians TCn and T^at. This 
allows one to expand the wave function of the A*"— reggeon state over the 
eigenstates of and as [91] 

^(6i, ...,bN) = J2Cab ^^"Hzi, zn)^^'Hzi, ...,zn), (4.23) 

a,b 

where Cab are the mixing coefficients. 
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Most remarkably, the Hamiltonian H4.2U() is in fact the Hamiltonian of the 
SL{2,C) Heisenberg magnet. The number of sites in the spin chain equals 
the number of reggeons and the corresponding spin operators are identified as 
six generators, L^,L^ and L^,L^, of the SL{2,C) group. The non-compact 
SL{2, C) Heisenberg magnet describes a completely integrable model. It 
possesses a large-enough set of mutually commuting conserved charges qn 
and Qn {n = 2,...,N) such that g„ = qli and [E[Ar,g„] = []H7v,^n] = 0. The 
charges Qn are polynomials of degree n in the holomorphic spin operators. 
They have the following form [7,8] 

Qn = ^jlj2^j2j3---^jnjlPjlPj2---Pjn (4.24) 

i<ii<j2<---<i„<iv 

with Zjk = Zj — Zk and pj = idz^ ■ The "lowest" charge q2 is related to the 
total spin of the system h. For the principal series of the SL{2,C) it takes 
the following values 

q2 = -h{h - 1) , h= + iuh , (4.25) 

with rifi integer and real. The eigenvalues of the integrals of motion, 
q2,...,qN, form the complete set of quantum numbers parameterizing the 
N—ieggeon states (|4.19|) . 

Identification of ()4.2U|) as the Hamiltonian of the 5L(2,C) Heisenberg 
magnet allows one to map the A^— reggeon states into the eigenstates of this 
lattice model. In spite of the fact that the noncompact SL{2, C) Heisenberg 
magnet represents a generalization of the compact SU{2) spin chain, very 
little was known about its energy spectrum until recently. The principal 
difficulty is that, in distinction with compact magnets, the quantum space 
of the SL{2, C) magnet does not possess a highest weight - the so-called 
"pseudo-vacuum state" - and, as a consequence, conventional methods like 
the Algebraic Bethe Ansatz method [32] are not applicable. The eigenprob- 
lem (|4.20() has been solved exactly in Refs. [27,28,92] using the method of the 
Baxter Q— operator [93] which does not rely on the existence of a= highest 
weight. In this approach, it becomes possible to establish the quantization 
conditions for the integrals of motion q^, ...,qisf and to obtain an explicit form 
for the dependence of the energy Ef^ on the integrals of motion. 

Thanks to complete integrability, the "conformal blocks" ^''•"^ and 
Eq. 1)4. 23() . diagonalize the integrals of motion in the holomorphic and an- 
tiholomorphic sectors, respectively. The function ^'(6i, . . . , b]\f) is a single- 
valued function on the two-dimensional plane. It belongs to the principal 
series of the SL{2, C) group, labelled by the spins (/i, h) defined in 1)4. 14() 
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(with h = 1 — h*), and is normahzable with respect to the SL{2,C) scalar 
product. In contrast with ^{bi, . . . ,bN), the chiral solutions ^^""^ and ^^^^ 
acquire nontrivial monodromy when encircles other particles on the plane. 
The quantization conditions for the integrals of motion q ensure that the 
monodromy cancels in the r.h.s. of (|4.2l-i|) . The same condition fixes (up to 
an overall normalization) the mixing coefficients Cab [28,91]. 

In this manner, the spectrum of the A^— reggeon state has been calculated 
for > 3 particles : For N = 3 few low-lying states have been found 
in [91,94] and the complete spectrum of states for 3 < iV < 8 was determined 
in [27,28] (see also [95]). Close examination revealed the following properties. 
The quantized values of the charges (with k = 3,...,N) depend on the 
"hidden" set of integers £ = {ii,i2, ■■■■,^2{N-2)) 

qk = qk{i^h;nh,(-) , (4.26) 

where the integer and the real number define the total SL{2, C) spin 
of the state, Eq. 1)4.14(1 . As a function of Vh, the charges form a family of 
trajectories in the moduli space q = {q2,Q3, ■■■,Qn) labelled by integers Uh 
and i . Each trajectory in the g— space induces a corresponding trajectory 
for the energy i?7v 

EN = EN{iyh;nh,£). (4.27) 

A few examples of such trajectories for = 3 are shown in Fig. |3 

Due to the complicated form of the quantization conditions, it remains 
unclear what the origin of trajectories is and what the physical interpreta- 
tion of the integers £ is. Both questions can be answered by solving the 
Schrodinger equation (|4.19|) within the semiclassical approach [96]. One 
might expect a priori that this approach could be applicable only for highly- 
excited states. Nevertheless, as was demonstrated in [96], the semi-classical 
formulae work with good accuracy throughout the whole spectrum. From 
the viewpoint of classical dynamics, the A^— reggeon state describes a chain 
of A^ interacting particles on the two-dimensional 6— plane [17,44,97]. The 
classical model inherits the complete integrability of the quantum noncom- 
pact spin magnet. Its Hamiltonian and the integrals of motion are obtained 
from (|4.2U|1 . (|4.18|) and (|4.24|) by replacing the momentum operators by the 
corresponding classical functions. Since the Hamiltonian ()4.2U|) is given by 
the sum of holomorphic and antiholomorphic functions, from point of view of 
classical dynamics the model describes two copies of one-dimensional systems 
"living" on the complex z— and z— lines. The solutions to the classical equa- 
tions of motion have a rich structure and turn out to be intrinsically related 



February 1, 2008 18:57 WSPC/Trim Size: 9.75in x 6. 5in for Proceedings BBGK 



324 A. Belitsky, V. Braun, A. Gorsky and G. Korchemsky 




-2-10 1 2 




Figure 5. The dependence of the energy E3 and the conserved charges qs on the total spin 
h = 1/2 + ii^Ti for rifi = 0. Three curves correspond to the trajectories with (<?i,^2) = (0, 2) , (2, 2) 
and (4,2). 

to the finite-gap solutions to the nonhnear equations [98,99]; namely, the 
classical trajectories have the form of plane waves propagating in the chain 
of A'^ particles. Their explicit form in terms of the Riemann 0— functions was 
established in [97] by the methods of finite-gap theory [98,99]. The charges 
q define the moduli of the finite-gap solutions and take arbitrary complex 
values in the classical model. Going over to the quantum model, one finds 
that charges q are quantized. In the semiclassical approach, their values 
satisfy the Bohr-Sommerfeld quantization conditions imposed on the orbits 
of classical motion of particles. 

In a standard manner, the WKB ansatz for the eigenfunction of the model 
(|4.2U() involves the "action" function, "^wKsi^i, ■ ■ ■ , zn) ~ exp(i5o/?i). Due 
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to complete integrability of the classical system, it can be defined clS 3) Sl~ 
multaneous solution of the system of the Hamilton-Jacobi equations 

f;g = P, q.(.,§^)=.„, (n = 2,...,iV), (4.28) 

where z = {zi, z^) denotes the set of holomorphic coordinates, qn{z,p) 
stands for the symbol of the operator ()4.24|) and P is a holomorphic com- 
ponent of the total momentum of N particles. The z— dependence of Sq is 
constrained by similar relations in the antiholomorphic sector. To find a 
general solution to Eq. (|4.28jl . one performs a canonical transformation to 
the classical separated coordinates [54, 98] 

(61,62, •••,&Af) {bo,xi,X2,...,XN-i) , (4.29) 

with 60 being the center-of-mass coordinate of the system, and x„ = 
{xn,Xn = x*) the new collective (separated) coordinates. Similarly to the 
example in Sect. 2.4, the classical dynamics in the separated variables is 
determined by the spectral curve ( "equal energy" condition) 

r^: = t%{x) - Ax^^ , (4.30) 

where y{x) = 2x^ sinhp^^ and px are the momenta in the separated coordi- 
nates. Here tAr(x) is a polynomial of degree N with the coefficients defined 
by the holomorphic integrals of motion g„ 

tN{x) = 2x^ + q2X^-'^ + ... + QN-lX + QN ■ (4.31) 

The spectral curve (|4.3U|) establishes the relation between holomorphic com- 
ponents of the separated coordinates, x and px, for a given set of the charges 

In the separated coordinates, the solution to the Hamilton-Jacobi equa- 
tions (|4.28|) takes the ioim So{zq,xi,X2, ■■■,xn-i) = {P ■ zq) + J2k=i ^oi^k) 
with [98] 

PX PX PX 

So{x) = / dxpx+ / dxpx = 2^e / dxpx- (4.32) 

J XO J XO J XQ 

Here px = Px is the complex momentum defined in (|4.3U|) and xq is arbitrary. 
The WKB expression for the wave function in the separated coordinates 
factorizes into a product of single-particle wave functions 

^(SoV) ^ ^^P■boQ^^^■^ . . QixN-i) (4.33) 



1,2008 18:57 WSPC/Trim Size: 9.75in x 6. 5in for Proceedings 



BBGK 



326 A. Belitsky, V. Braun, A. Gorsky and G. Korchemsky 

with Q{xk) ~ exp (z5o(xfc)). According to (|4.3U|) . the momentum, px, and, 
as a consequence, the action function So{x) are multi- valued functions of x. 
Denoting the different branches of the action function by 5o,cj(x), one writes 
the WKB expression for the wave function of the quantum spin magnet as 
a sum over branches [17, 100] 

Q{x) = Y,Mx)exp Uso,a{x)] , (4.34) 

a ^ ^ 

with h = 1. The function Ak{x) takes into account subleading WKB correc- 
tions and is fixed uniquely by So^aix)- The quantization conditions for the 
charges q follow from the requirement that (|4.34jl has to be a single-valued 
function of x. As was shown in Refs. [44,96], these conditions can be ex- 
pressed in terms of the periods of the "action" differential over the canonical 
set of the a— and /9— cycles on the Riemann surface corresponding to the 
complex curve ()4.30l) 



d> dx px = TT £2k-i , 'Sie d> dxpx = TTi2k, (4.35) 

with k = 1, .., N — 2 and £ = (£i, . . . , l2N-i) being the set of integers. The 
relations 1)4. 35(1 define the system of 2(A^ — 2) real equations on the (A^ — 2) 
complex charges q3,...,qN (we recall that the eigenvalues of the "lowest" 
charge q2 are given by (|4.25|) ). Their solution leads to the semiclassical 
expression for the eigenvalues of the conserved charges. In turn, the energy 
of the A'— reggeon states -Ejv,q can be expressed as a function of qs, ...,(?Ar. 
In the semiclassical approach, the corresponding expression is 

^5^'^ =4 In 2 (4.36) 



N r 

+ 2Ke^ 

fc=0 



-0(1 + meXk + l^mAfcl) + V(^KeAfc + |9mAfc|) - 2V^(1) 



where are roots of the polynomial iAr(u) (see Eq. ()4.31|) '). The expression 
in Eq. (|4.36|1 is similar to the energy of the SL{2,M) magnet in Eq. (|2.69j) 
although the properties of the two models are different. As was demonstrated 
in [96], the resulting semiclassical expressions for q3,...,qN and Ejy are in 
good agreement with exact results [27,28]. As an example, we show in Fig. 
a comparison of the exact spectrum of the conserved charge q^ for the A^ = 3 
states carrying the conformal spin h = 1/2 (or equivalently = i^h = 0) 
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Figure 6. Lattice structure at Af = 3. Crosses denote the exact values of q^' at h = 1/2. Dotted 
lines intersect at the points defined in Eq. 14.371 . 



with the semiclassical expression given by 

r3(2/3) 



1/3 
93 



2tt 



(4.37) 



with integer li and £2 introduced in (|4.26|) . 

A novel feature of the quantization conditions (|4.35j) is that they in- 
volve both the a— and /3— periods on the Riemann surface. This should be 
compared with the situation in the SL(2,M) Heisenberg magnet discussed 
in Section 2.4. There, the WKB quantization conditions involve only the 
a— cycles, Eq. ()2.67|) . since the /?— cycles correspond to classically forbidden 
zones. For the SL{2, C) magnet, the classical trajectories wrap over an ar- 
bitrary closed contour on the spectral curve H4.3U() leading to H4.35() . This 
fact allows one to explore the full modular group [101] of the complex curve 
H4.3U() and explain the lattice structure of the spectrum shown in Fig. El 
We refer the interested reader to Ref. [96] for details. 



5. Spin chains and Af = 2 super Yang-Mill theories 

In this section, we shall briefly discuss the integrability phenomenon which 
emerged from the calculation of the low-energy effective action in M = 2 
SYM theory in Refs. [29, 102]. Rather than delivering an extensive review of 
the subject, we shall outline the similarities and differences with the cases 
discussed in the previous sections. Detailed reviews of the solutions toM = 2 
SYM theories are given in [103, 104] while reviews of the integrable struc- 
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tures in low-energy effective actions can be found in [105-107] . We will argue 
below that the low-energy effective action in the M = 2 theory with matter 
in the fundamental representation is governed by integrable models which 
belong to the same universality class as the spin chains that we encountered 
in the studies of anomalous dimensions of Wilson operators and high-energy 
scattering in QCD. The particular form of the model depends on the minibcr 
of matter hypcrmultiplcts Nf. To preserve the asymptotic freedom in the 
A/" = 2 SYM theory with the SU {N^) gauge group, one requires that Nj does 
not exceed 2Nc. For Nf = 2Nc the beta-function in the M = 2 theory van- 
ishes and the low-energy effective action is described by an inhomogeneous 
twisted Heisenbcrg XXX spin chain while for Nf < 2Nc one encounters its 
degenerations [108,109]. 

The classical action of the J\f = 2 SYM theory has the following form 



where r is the complex coupling constant 



$ is the scalar, chiral M = I superfield involving the scalar fields and 
fermions. V is the vector = 1 superfield involving the gauge field and 
its superpartner (the gluino), W is the spinor superfield constructed from 
V in such a way that trp^^ gives rise to kinetic terms for the gluons and 
gluinos. The superfields Qk and Qk (with k = 1,. . . ,Nf) represent matter 
in the fundamental and anti-fundamental representations of the gauge group 
SU{Nc), respectively, while the superfields $ and W belong to the adjoint 
representation. The J\f = 1 superfields W and <1> can be combined together 
into the M = 2 chiral superfield '^j\f=2, while pairing Qk and Qk yields the 
M = 2 hyper multiplet. 

In the N = 2 SYM theory, the vacuum states which do not break super- 
symmetry arc field configurations with vanishing energy. It turns out that 
there are vacuum valleys, parameterized by the vacuum expectation values 
of the fields Qk and We shall be interested in the so-called Coulomb 
branch where all matter condensates vanish, [Qk) = 0, and the equation 
defining the condensate of the scalars reads 




(5.1) 



tr[</,,</,t]2 = 0. 



(5.3) 
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Up to an SU{Nc) gauge rotation, the expectation values of the scalar field can 
be chosen to lie in the Cartan subalgebra 4> = diag(ai, ...,aN^) with ti(p = 
0. The parameters cannot serve, however, as good order parameters, 
since there is still a residual Weyl symmetry which changes but leaves 
gauge-invariant quantities intact. An appropriate set of gauge invariant 
order parameters that fix the vacuum state unambiguously is provided by 
Uk = (tr (j)'') with k = 2, . . . , Nc- The variables U2, ■ ■ ■ ,un^ parameterize the 
space of vacua in the M = 2 theory - the moduli space, and the choice of a 
point on the Coulomb branch is equivalent to a choice of the vacuum state. 
Simultaneously, the expectation value of the scalars yields the scale at which 
the coupling constant is frozen. At a generic point on the moduli space, 
for non-vanishing (j), the SU{Nc) gauge symmetry is broken by the standard 
Higgs mechanism and hence massive gauge bosons emerge. Integrating them 
out, one arrives at the well-defined problem of the determination of the low- 
energy effective action of the N = 2 SYM theory. 

To calculate the M = 2 low-energy effective action, one has to modify 
the classical action (|5.1() by taking into account one-loop perturbative cor- 
rections and, in addition, include instanton corrections. The explicit sum- 
mation over the instantons is a complicated problem which has been solved 
only recently [110]. In their original work [29, 102], Seiberg and Witten 
avoided this problem and calculated the low-energy effective action by mak- 
ing use of three ideas: holomorphicity, duality and their compatibility with 
renormalization-group flows. Holomorphicity implies that the low-energy ef- 
fective action depends on a single holomorphic function T called prepotential 



Thus, the problem effectively reduces to the determination of the prepo- 
tential. The holomorphic function J-{a) is defined unambiguously by its 
behavior in the vicinity of singular points on the complex a— plane. One of 
the singularities, corresponding to large values of the condensate, is fixed 
by the perturbative contribution to the M = 2 low-energy effective action. 
The behavior at the remaining singularities is fixed by the use of duality and 
non-renormalization theorems for the central charges of the SUSY algebra. 

To define the duality transformations relating the weak and strong cou- 
pling constant regimes in the M = 2 SYM theory, it is instructive to examine 
first the M = A model. In this model, the coupling constant does not run 
and the duality transformations are generated by r — > and r — > r + 1. 
The same transformations define the modular group acting on the moduli 




(5.4) 
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space of some Riemann surface parameterized by r. This suggests that the 
duahty properties of the SYM theory are encoded in the properties of Rie- 
mann surfaces. Going over from M = 4 to M = 2 theory, one finds that 
naive generahzation of the duahty meets serious difficulties. The reason for 
this is that, in the asymptotically free theory, one has to marry the duality 
with the renormalization group flow of the coupling constant r. The solution 
to this problem is that one is still able to connect the duality and modular 
transformations if one considers the M = 2 theory at different vacua, con- 
nected to each other by the duality transformation. Then, the duality acts 
on the moduli space of vacua which is associated with the moduli space of 
the auxiliary Riemann surface. At the next step, one has to find out proper 
variables whose modular properties fit the field theory interpretation. These 
variables are identified with the integrals of a meromorphic 1-form dS over 
cycles on some Riemann surface [29, 102] 



where i = 1, Nc — I for the SU (Nc) gauge group. 

The integrals 1)5. 5() play a two-fold role in the Seiberg-Witten approach. 
Firstly, they allow one to determine the prepotential J- (a) through the rela- 
tion 



with a = (ai, . . . , aj\fc-i) and, therefore, evaluate the low energy effective ac- 
tion 1)5.4(1 . Then, using the property of the differential dS that its variations 
w.r.t. moduli are holomorphic, one can also calculate the matrix of coupling 
constants 



Secondly, the periods o and ao, Eq. (|5.5|) . define the spectrum of stable states 
in the theory which saturate the Bogomolny-Prasad-Sommerfeld (BPS) 
bound. For instance, the dyonic spectrum of stable particles of a SU{2) 
theory reads 




(5.5) 




(5.6) 




(5.7) 





with n and m integeral. 

The Seiberg-Witten solution admits an elegant interpretation in terms 
of classical integrable systems. The Riemann surface entering H5.5|) defines 
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the solution to the classical equations of motion in a classical integrable 
system. Moreover, the meromorphic differential dS turns out to be the action 
differential in the separated variables in the same system, dS = p{x) dx, 
with p{x) being a single- valued function on the Riemann surface [111]. The 
dependence of the prepotential on the fundamental scale A is given by [112] 

= Pu2ia) ^ (3H , (5.9) 

where /3 is the one-loop beta-function in the N = 2 SYM theory. This equa- 
tion has another interpretation as the evolution equation in an integrable 
system. H coincides with the Hamiltonian of this system and In A is the 
evolution time variable. Hence we see, once again, that the logarithm of the 
relevant scale plays the role of time in the integrable dynamics. 

As was already mentioned, the Riemann surface parameterizing the 
Seiberg-Witten solution for the low-energy effective action in the M = 2 
SYM theory with fundamental matter [108] coincides with the spectral curve 
describing the solution to the classical equation of motion in the periodic in- 
homogeneous Heisenberg XXX spin chain of length [113,114]. The masses 
of the squarks in the N = 2 theory are mapped into the parameters of the 
chain while the coordinates on the Coulomb branch of the moduli space are 
mapped into the values of the integrals of motion. The spectral curve for 
the spin chain system is defined as 

det('u; -rAr^(x)) = 0, (5.10) 

where T]\f^{x) = tr[Li(3;) . . . Ljv^(x)] is the transfer matrix of the model and 
Li{x) is a local 2x2 Lax matrix given by 



L,(x) = (x + Ai)•l + ^5a,^•a^ (5.11) 



a=l 



Here cr" are the Pauli matrices and Aj are the chain inhomogeneities. The 
Poisson brackets of the dynamical variables S'a, (a = 1, 2, 3) are just 

{Sa,i,Si,^k} = —i^abcScjSjk , (5-12) 

so that the vector {Sa} plays the role of the angular momentum ("classical 
spin"). The spectral curve (|5.10j) is now 

w + ^^ = 2P{x), (5.13) 
w 
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with P{x) = tiT{x)/2 and Q{x) = detr(x). In the hyper-eUiptic parame- 
terization, y = {w + Q{x) /w) /2, the same curve becomes 

The zeroes of Q{x) define the masses of the hypermultiplets in the M = 2 
SYM theory via 

Q(x) = JJdetLi(x) =Yl{x-mt){x-m-) (5.15) 

1=1 i=l 

where mf = —Xi±Si, = sf and Si stands for the spins at the feth 
site of the chain. By construction, the number of matter multiplets is twice 
the number of colors. Still, Eq. 1)5. 14() is not the most general form of the 
spectral curve for Nf = 2Nc. 

Equation (|5.14|) can be generalized by taking into account the additional 
freedom in the definition of the spin chain and the spectral curve. Namely, 
one can multiply the Lax operator of the spin chain by an arbitrary constant 
matrix without changing the commutation relations and conservation laws. 
Moreover, one can also attach a constant (external magnetic field) matrix 
V to the end of the chain, that is L]y^{x) Ln^{x)V. This amounts to 
imposing more general (twisted) boundary conditions along the spin chain. 
The spectral curve corresponding to the "twisted" spin chain takes the form 
[102,113,114] 

w + ^^ = P{x), (5.16) 




2/i \ 

-rrij > , "T-j • 

Let us now consider these formulae in the special case of a superconformal 
N=2 SYM theory with Nf = 2Nc massless fundamental hypermultiplets 
[102, 113-117]. The corresponding integrable system is described by the 
spectral curve 

y^ = P\x)-Ax^''^{l-p\T,i)), (5.17) 

where {t^i ) is some function of the coupling constant in the M = 2 theory 
and the polynomial P{x) depends on the coordinates on the moduli space 
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U= {U2, ■■■,UnJ 

= J2 Iki^) ^"^^"^ = 2^"^^ + 1^ ^"^^"^ + - + IN. ' (5-18) 

A:=0 

where Qq = 2, qi = and the other qi. are some known functions of u. 
Their exphcit form is not important for our purposes. The Seiberg-Witten 
meromorphic differential on the curve is given by 

dS = pdx = ln{uj/x^'')dx, (5.19) 

where y = oo — x'^^" /to and to = x^"^ . 

Prom the point of view of integrable models, the spectral curve 1)5. 17() 
corresponds to a classical Heisenberg XXX spin chain of length Nc with 
spin-zero at all sites (due to qi = 0) and parameter p related to the exter- 
nal magnetic field, or equivalently, to the twisted boundary conditions [109]. 
Remarkably enough, the spectral curve 1)5. 17(1 describing the low-energy ef- 
fective action of the superconformal N=2 SYM theory with Nf = 2Nc is 
closely related to the spectral curve ()4.3U() describing the spectrum of multi- 
reggeon compound states in multicolor QCD [44]. It is easy to verify that 
the two curves coincide if one makes the following identification: 

• The number of the reggeons N = N^, 

• The integrals of motion for the multi-reggeon state are identified as 
the above-mentioned functions qkiu) on the moduli space of the su- 
perconformal theory; 

• The coupling constant of the gauge theory should be such that p{tci) = 
that is Tci = ^ + I tan 

Under these conditions the two theories fall into the same universality class. 

Let us briefly mention some recent developments concerning the relation 
between integrable systems and M = 1 SYM theories. In the N = 2 SYM 
theory, the Riemann surface degenerates at some points on the Coulomb 
branch of the moduli space. After soft breaking oi N = 2 down to N = 1, 
these points correspond to the vacuum states of the N = 1 SYM theory. 
Some massless states condense at these points leading to the formation of a 
mass gap and to confinement. For instance, in the SU{2) case the M = 1 
vacua correspond to the points U2 = =bA^ where the monopoles or dyons be- 
come massless and condense. It turns out that the relation to the integrable 
systems becomes even more direct after soft breaking. Namely the tree- 
level superpotential Wtree = X^fc tj^tr^'' amounts to the following effective 
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superpotential 



(5.20) 



k 



where L is the Lax operator of the corresponding integrable system gov- 
erning the M = 2 theory with the same matter content [118 120]. That is 
^eff = '^k^kHk, and vacua in tlic J\f = 1 theory defined by the extrem- 
ization of Wes are in one-to-one correspondence with the equilibrium states 
of the corresponding integrable system with respect to the combination of 
Hamiltonians above. 

6. Gauge/string correspondence 

In previous sections, we described the hidden integrability symmetry of the 

one-loop dilatation operator in QCD and its supersymmetric extensions. Its 
origin, however, remains obscure on the gauge theory side. In this section 
we shall describe attempts to explain the integrability phenomena in Yang- 
Mills theories based on the gauge/string correspondence. We will demon- 
strate that quantum integrability of the dilatation operator in gauge theory 
is ultimately related to classical integrability of the relevant stringy sigma 
models. The subject is rapidly developing but a fully satisfactory picture 
has not yet emerged. 

To start with, let us recall the general features of the gauge/string cor- 
respondence relevant for our discussion. In the case of the Af = 4 SYM 
theory it states that this gauge theory is dual to IIB string theory on an 
AdS^ X background with additional flux of a higher-form field [38, 121]. 
The background metric in Poincare coordinates has the form 



and can be considered as the near-horizon limit of the D3 brane metric. 
Since the D3-brane is a source of the Ramond-Ramond four- form field A4, 
the background solution is supplemented by the flux of the corresponding 
field strength 



ds^ = ^{-dt^ + dxj + dxl + dxl) + R^-^ + B?d^l 



(6.1) 





(6.2) 



The radii of the AdS^ and are identical and equal to 



= 4.Trgsa'^Nc . 



(6.3) 
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The four-dimensional gauge theory is locahzed on the boundary of the AdS^ . 
The conformal 5*0(2, 4) group and the _R— symmetry 50(6) group of the 
M = ^ SYM theory are identified with the isometry group of the AdS^ and 
spaces, respectively. According to the gauge/string duality, the eigenvalues 
of the dilatation operator in the A/" = 4 SYM theory and the energy spectrum 
of the string in radial quantization coincide. 

The string tension is related to the gauge theory 't Hooft coupling con- 
stant by 

Therefore the strong coupling regime A ^ 1 in the gauge theory corresponds 
to the semiclassical regime on the stringy side. In this way, the classical string 
calculations provide predictions for many interesting quantities in the A/" = 4 
SYM theory (like anomalous dimensions, condensates, Wilson loops etc.) in 
the strong coupling regime [121]. However, these predictions can be checked 
diagrammatically only in a few examples, such as for a circular Wilson loop 
[122, 123]. On the other hand, to approach the weak coupling regime in 
Yang-Mills theory, one should go beyond the semiclassical approximation 
and take into account quantum effects on the string theory side. 

At present, an explicit quantum solution of the string theory in the 
AdS^ X background is not available. This makes it impossible to compare 
the string spectrum with the complete set of operators in A/" = 4 SYM the- 
ory. Hopefully the explicit solution of the string theory can be found [124] 
in the case of the limiting geometry defined as the Penrose limit of the 
AdS^ X [39,125,126]. The Penrose limit can be described as the region 
around the null geodesic in the AdS^ x . Introducing new variables 

x+ = ^, x-=iiR\t-x) (6.5) 

with X being an angular variable in and /i being some scale, one takes 
the limit R—*co and recovers the pp-wave metric 

8 

ds"^ = -Adx+dx- - z^dx+^ + ^ dzf . (6.6) 



Here, eight, flat, transverse coordinates Zi come both from the and AdSr, 
parts of the geometry. In this metric, the string behaves as a particle rotating 
with large angular momentum J along the angular coordinate % in 5^. The 
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light-cone energy of the string is 

H = 2p- = i{dt + d^) = (A - J) . (6.7) 

For i? — > oo it takes finite values provided that the following double scaling 
limit is considered: 

R ^ oo , A~J^oo, = const . (6.8) 

Quantization of the string propagating in this background reduces to the 
quantization of the oscillators. As a result, the exact spectrum of the type 
IIB string in the pp-wave background looks like 

^-'^ = E^^\A + ^' (6.9) 

where k labels the Fourier modes, Nj^. denotes the total occupation number 
of oscillatory mode and the condition P = kN/. = is imposed. 

On the gauge theory side, (|6.9j) defines the anomalous dimensions of 
certain Wilson operators in the = 4 SYM theory. The length of the 
string J equals the number of constituents of the composite operator. The 
ground state of the string can be identified with the operator built from 
scalars Z = $i + 

\0,J)^trZ-^. (6.10) 

These operators have charge J with respect to the "rotation plane" in the 
pp-wave. The oscillatory excitations of the string ground state correspond to 
insertions of other scalar fields. For example, the so-called BMN operators 
in the = 4 SYM theory can be mapped into the stringy modes as follows 

4+|0, J) -^tr^>iZ'^, 

a;+ai+ JO, J)^Y1 e^^'/'^tr ^iZ'^^Z-^"' . (6.11) 
I 

One can deduce from these expressions that calculation of the stringy spec- 
trum corresponds to diagonalization of the mixing matrix for Wilson oper- 
ators on the gauge theory side. 

The energy of the string in the pp-wave limit is a function of the ratio 
of the coupling constant and angular momentum, A/J^. It is expected that 
expansion of this function in powers of A/ should reproduce a perturbative 
series for the anomalous dimension of the corresponding Wilson operators in 
the weak-coupling regime. On the other hand, the one-loop J\f = 4 dilatation 
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operator in the sector of scalar operators coincides with the Hamiltonian 
of the Heisenberg SO{6) spin chain [25]. This allows one to map stringy 
states into spin chain states. The correspondence is very precise for Wilson 
operators built only from two complex scalars, in which case the iS'0(6) spin 
chain reduces to the conventional Heisenberg SU{2) spin— 1/2 chain. Then, 
the ground state in the string theory corresponds to all spins aligned in the 
same direction in the isotopic space while the stringy excitations correspond 
to flipping some spins along the chain. 

The BMN operators are the only well-established example in which the 
exact quantum answer on the stringy side can be matched into the all-loop 
anomalous dimension on the gauge theory side. Later in this section, we 
shall discuss generalized BMN operators which can be treated semiclassi- 
cally on the stringy side and demonstrate that the corresponding solutions 
to the classical equations of motion in the stringy sigma models are ulti- 
mately related to classical integrable models. We shall explain the relation 
between the string sigma model and quantum spin chains in the thcrmody- 
namical limit and demonstrate that the anomalous dimensions of the certain 
operators are in one-to-one correspondence with the special class of classical 
solutions to the sigma model for which this model reduces to finite dimen- 
sional integrable systems of the Neumann type. Finally, we shall comment 
on the relation between general classical solutions to the sigma model and 
the Bethe Ansatz solution to the compact quantum spin chains in the semi- 
classical limit. 

6.1. Derivation of the string in the thermodynamical limit 

The sigma model describing the string moving in the appropriate curved 
background can be derived from the quantum compact spin chain in the 
long- wavelength limit. The corrections to the classical sigma model scale 
as 1/J, where the angular momentum of the string J corresponds on the 
gauge theory side to the number of fields entering the composite operator 
or equivalently the length of the spin chain. The transition from the spin 
chains to the sigma model relies on the coherent state formalism [42] . 

Let \ss) be the state with total spin s and projection onto the z— axis 
Sz = s. The coherent state for the spin— s representation of the SU{2) 
group is defined as 

\fi) = e'^-^ e'^y^ \ss) (6.12) 
where n is the unit vector, = 1, 

n = (sin 9 cos 4>, sin 6 sin 0, cos 6) (6.13) 
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with 6 and (p being spherical angles. Expanding the Hamiltonian of the 
spin chain H = X/ (47r^) X]fc=i(l/4 — Sk ■ Sk+i) over the coherent states, one 
rewrites the partition function tre~^* in the standard manner as a path 
integral over = sftk with the following action 

5(n) = I dt I dTftkidtfik x drUk) - ^■^^ / dt'^{nk - fik+if , 

k=l J J 

(6.14) 

with nj_|_i = ni. In the long- wavelength limit, the vectors nk{t) vary 
smoothly along the spin chain and, therefore, they can be approximated 
by a function n(c7, t) with continuous a running between and the chain 
length J, leading to 

dtda dt(j) cos e - -AtS^ I dtda \idjf + id^^f O] . (6.15) 
Stt^ J 

It turns out [42] that for s = 1/2, this expression coincides with the stringy 
action 

5str = ^ / <i^dT [G^^drX'^drX'^ - G ^.O^X'^d^X'^] (6.16) 



evaluated for the classical string propagating in the background ds 



2 



Gf,^dX^'dX' 



ds'^ = -dt^ + dil? + difl + dipl + 2 cos{2ilj)dipidip2 . (6.17) 



To see this, one fixes the gauge t = x^i takes the limit drX^ and x ^ oo 
with xdrX^ = fixed and identifies the variables as 

(^2 = -^<^, ^ = 1^- (6-18) 

Then one eliminates (pi with the help of the classical equations of motion 
and arrives at (|6.15j) . 

The derivation of the effective action can be also generalized to the SU (3) 
case [127] and to the string carrying both large Lorentz spin S and i?— charge 
J [128]. One can improve the effective sigma model action (j6.15|) by calcu- 
lating corrections involving higher derivatives of the fields. All such terms 
containing up to four derivatives have been found in Refs. [129,130]. 



6.2. Semiclassical string motion and integrable models 

We have argued above that the effective action for long-wavelength exci- 
tations in the compact spin chain coincides with the classical action of the 
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sigma model on the curved background relevant for calculation of the anoma- 
lous dimensions of BMN-like operators. As the next step, the corresponding 
solutions to the equations of motion are compared. To this end, one consid- 
ers the bosonic part of the superstring action on the AdS^ x background. 
It is given by the sum of two coset sigma models 

^ = ^1 dadrlGifdyA + G^ax'^dx'] , (6.19) 

where the string tension is proportional to the 't Hooft coupling, Eq. (|6.4j) . 
It is convenient to rewrite the action with the constraint imposed by the 
Lagrangian multiplier 

S = ^ j dadT[dXmdXm + Ax{X^ - 1) + dY'^dV'' + Ay{Y^ + 1)] (6.20) 

where Xn {n = 1, ... 6) and Yk {k = 0, . . . , 5) are the two sets of the em- 
bedded coordinates in the flat space with signatures (6,0) and (4,2), 
respectively. The action has to be supplemented by the Virasoro constraint 
for the vanishing of the two-dimensional energy momentum tensor 

YkYi + YUl + XnXn + X;x; = Yj,Yi + X^X'n = (6.21) 

and by the periodic boundary conditions 

Yk{a + 2^)=Yk{a), X„(a + 2^) = X„(a) . (6.22) 

Due to the 50(2,4) and 5*0(6) symmetries, the classical action possesses 
the set of conserved charges 

Ski = VxJ da{YkYi - YiYk) , 

Jnm = J da{XnXm — X^Xi) . (6.23) 

Among them, one distinguishes 6 Cartan generators: the energy E = Sq^, 
the Lorentz spins 5i2, 5*34 and the angular momenta J12, J34, Jse- These 
conserved charges parameterize general solutions to the classical equations 
of motion [131]. 

To describe a particular operator on the gauge theory side, we have to 
identify the corresponding solution to the classical equations of motion in 
the sigma model (|6.20p . subject to the constraints (|6.21jl and (|6.22jl . The 
simplest ansatz, corresponding to a string located at the center of the AdS^ 
and rotating in the S^, is 

Y, + iYo = e^* , X2i-i + iX2i = r^ia) e'^-^^+'^^i-) . (6.24) 
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with i = 1,2,3 and the remaining 1"— coordinates set to zero. Substituting 
into the sigma model action (|6.2U() yields the Lagrangian [132] 

^ = E {r? + rfc^ - -^M) - A. ^ (rf - l) . (6.25) 

i=l i=l 

Solving the equations of motion for Oi one gets a[ = Vi/rf with Vi being 
the integration constants. The resulting Lagrangian describes the integrable 
Neumann-Rosochatius system.' It admits five independent integrals of mo- 
tion: ^1,^2,^3 plus two additional integrals of the form 



3 ^ 

li ='^i ,,,2 _ ,.,2 



(6.26) 



subject to Yl\=i li = 0. The periodicity condition trades Vi for three inte- 
gers rrii and li for two integers nj. As a result, the energy depends on the 
frequencies and five integers. These variables are not independent since 
the Virasoro constraint imposes a relation between them. For this type of 
string motion, the infinite set of conserved charges in the sigma model is pa- 
rameterized by a finite set of integrals of the motion [133, 134]. The energy 
corresponding to classical solutions of the string sigma model defines the 
anomalous dimension of the dual composite scalar operators in the = 4 
SYM theory. There are many examples of such correspondences discussed in 
the literature, initiated in Ref. [131] and further developed in Refs. [135-140]. 

It is interesting to note that the Neumann system is isomorphic to the 
stationary solutions of the Landau-Lifshitz equation. We recall that the 
time variable in the Neumann system is identified with the coordinate along 
the string. This identification survives time discretization in which case the 
Neumann system is equivalent to the discretized version of the Landau- 
Lifshitz equation describing the XYZ spin chain [141] and it can be used for 
the formulation of the integrability in the context of the string bit model 
[142]. 

In the above example, the string was rotating in the S^. In general, it 
moves both in the AdS^ and and could have large angular momenta in 
both spaces. Contrary to the situation with pure scalar operators, when the 
comparison with the loop expansion on the gauge side can be performed for 
operators with large R charge, the situation with operators carrying large 
Lorentz spin S is more subtle. The folded, closed string rotating in the AdS^ 



' For vanishing Vi the system 16.251 reduces to the Neumann model with three degrees of freedom. 
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yields the dependence for the anomalous dimensions of twist-two operators 
F-^-±{D^)^ F^± at large coupling of the form [40] 

^^-=^) = :0lnS\ (6.27) 

This result can be generalized to higher twist operators of the form 
F^±D_^} Fj^± . . . D_t_^~^ F^±. The energy of the corresponding revolving 
string coincides with the energy of the classical Heisenberg spin chain of 
length L and leads to [143] 

7t£U-. = ||lngL(Si,52,...,5^_i). (6.28) 

Here qi is the highest integral of motion of the spin chain. For S'fc ~ 5" S> 1 
with k = 1, . . . , L — 1 one has ql ~ S^. 

Notice that the logarithmic scaling of the anomalous dimensions is a uni- 
versal feature of Wilson operators with large Lorentz spin in gauge theories, 
unrelated to the presence of supersymmetry [144,145]. However stringy de- 
scription of this scaling at weak coupling remains unknown and it is doubtful 
whether such classical string sigma model solutions exist [140]. One should 
expect instead that ~ In behavior at weak coupling is driven by the quan- 
tum sigma model. 

The integrability phenomenon offers the possibility of extending the 
gauge/string duality beyond the special class of classical string solutions 
described above. Namely, instead of comparing particular solutions one can 
identify the integrable structures corresponding to quantum spin chains de- 
scribing the dilatation operator on the gauge theory side and to classical 
equations of motion on the string theory side [146]. It turns out that, in 
both cases, integrability is encoded in the properties of Riemann surfaces. 

For the quantum spin chains, the appearance of Riemann surfaces within 
the framework of the Bethe Ansatz is not surprising. As we already explained 
in Sect. 2.4, semiclassical solutions to the Baxter equation are determined 
by the properties of the spectral curve whose genus is proportional to the 
number of sites in the chain. In Sect. 2 we discussed baryon operators built 
from three quarks and the corresponding Riemann surface, Eq. (|2.66|1 . had 
genus equal to 1. In the case of the BMN-like operators, the number of 
constituent scalar fields goes to infinity in the thermodynamic limit and, 
therefore, the corresponding Riemann surface would have an infinite genus. 
However, by choosing the appropriate values of the integrals of motion, it 
is possible to pinch almost all handles and obtain a finite genus surface. 
It is this degenerate surface which parameterizes general solutions to the 
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classical equations of motion of the string [146,147]. The agreement between 
semiclassical solutions to the Bethe Ansatz equations and solutions to the 
classical string equations of motion has been carefully checked up to the 
two- loop level [146]. 

Although the correspondence between one- and two-loop ™ dilatation 
operators in the = 4 SYM theory Yang-Mills theory, stringy states 
and integrable quantum spin chains is well established, the situation with 
higher loops in perturbation theory is unclear. Several proposals have been 
made concerning integrable structures behind a higher loop dilatation opera- 
tor [148,149]. At the same time, starting at three-loop order the discrepancy 
seems to arise between expressions for the anomalous dimensions of compos- 
ite scalar operators with large-i? charge and the energy spectrum of the 
string [150]. More work is needed to clarify this issue. 



6.3. Open string picture for anomalous dimensions 

There exists an alternative description of the logarithmic growth of 
the anomalous dimensions of Wilson operators with large Lorentz spin, 
Eq. (|6.27|) and (|6.28|) . in terms of Wilson loops in the gauge theory and 
open strings on the AdS^ background. This picture relies on the correspon- 
dence between the anomalous dimensions of the composite operators with 
large number of light-cone derivatives and the so-called cusp anomaly of 
Wilson loops [144, 145]. It was shown a long time ago [151] that the Wil- 
son loop W[C] = tr{Pexp(i(7 dx'^^^(x))} acquires a nontrivial anomalous 
dimension TcuspC^, 0) if the integration contour C has a cusp 

{W[C]) ~ /-=p(^'^) , (6.29) 

with fi being a UV cut-off. The cusp angle 6 is restricted to the interval 
[0, 27r[ in Euclidean space but is unrestricted in Minkowski space. The corre- 
spondence between the anomalous dimension of twist-2 spin operators with 
large Lorentz spin S and the cusp anomaly is follows [144, 145] 

^^tw=2) ^ 2reusp(A, e = In S) (6.30) 

and is valid for an arbitrary coupling constant A. At weak coupling and 
9^1, one has 



rcusp(A, 6) 



(6.31) 



In a closed subsector only. 
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with perturbative coefficients known up to three-loop order [152]. The cal- 
culation of rcusp(A, 6*) at the strong coupling can be effectively done via the 
open string picture. In this limit, the Wilson loop is proportional to the area 
of the minimal surface swept out by an open string which penetrates into 
the fifth AdS dimension and whose ends trace the integration contour C in 
Minkowski space [122]. This leads to [153,154] 



for ^ ^ 1. Being combined together, Eqs. (|6.3U|) and (|6.32|) reproduce the 
strong coupling result ()6.27|) based on the folded closed string picture [40] . 

The correspondence 1)6. 3U() can be generalized to higher twist operators. 
In that case, the anomalous dimension of the Wilson operator built from L 
constituent fields and a total number of derivatives S, such that S ^ L, can 
be mapped into the anomalous dimension of the product of L Wilson loops 
in the fundamental representation of the SU{Nc) group; the total number of 
cusps varies between 4 and 2L [143]. At large Nc, the expectation value of 
the product of Wilson loops factorizes into the product of their expectation 
values. This implies that, at strong coupling the area of the minimal surface 
corresponding to the product of /c = 2, . . . , L Wilson loops with cusps is 
given by the sum of k elementary areas leading to 



2rcusp(A,0 = InS) < 7r^'^^(A) < Lr,^,p{X,9 = \nS) . (6.33) 



We are reminded that the anomalous dimensions of higher twist operators 
are not solely determined by the total number of derivatives S. They form 
instead a band whose internal structure at weak coupling is governed by 
integrals of the motion of the quantum Heisenberg SL{2) magnet. Equation 
(|6.33|) defines the boundaries of the band both at strong and weak coupling. 

7. Conclusion 

The main objective of the present work was to demonstrate how the phe- 
nomenon of integrability arises in certain limiting cases of Yang-Mills dy- 
namics and exhibit their similarities or ultimate relation whenever it was 
obvious. In all cases we have studied, integrability appears as a hidden sym- 
metry of an underlying effective theory. It either describes elementary fields 
living on the light cone, in the case of renormalization group evolution in 
QCD and its supersymmetric extensions, or reggeons as new degrees of free- 
dom in high-energy QCD, or branes as effective degrees of freedom for the 
Seiberg-Witten solution to the M = 2 SYM theory. A number of nontrivial 




(6.32) 
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questions still have to be answered. The most obvious and at the same time, 
the most profound and most difficult question is "What is the origin of inte- 
grability?" or, in other words, "What is the symmetry, if any, of the gauge 
theory which leads to it?" More specific questions that have recently kept 
theoretical physicists busy concern issues like the fate of the integrability 
of the dilatation operator in SYM at higher orders of perturbation theory, 
the integrability of the SYM dual sigma models on curved backgrounds and 
the matching of integrable structures on both sides of the correspondence, 
to name but a few. Ultimately, if the integrability is indeed a property of 
the full quantum Yang-Mills theory as well as the dual string theory, it will 
provide the most sophisticated test of duality and endow us with powerful 
machinery to tackle the strong-coupling regime of field theories. On this 
quest, unfortunately, we will be missing Ian Kogan. 
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